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ABSTRACT 

U.  S.  Navy  ships  not  equipped  with  NTDS  currently  per¬ 
form  surface  radar  target  tracking  manually,  a  tedious  and 
inaccurate  process. 

This  thesis  establishes  that  an  8-bit  microprocessor, 
interfaced  to  a  common  surface  search  radar,  could  compute 
the  course  and  speed  of  multiple  radar  targets.  The 
investigation  includes  the  estimation  of  radar  measurement 
errors,  a  determination  of  practical  digital  interface 
limitations,  and  the  development  of  a  Kalman  filter 
algorithm,  a  floating  point  arithmetic  library  for  the 
INTEL  8080  microprocessor,  aind  the  microprocessor  tracking 
system  software.  The  error  statistics  and  execution  time 
of  the  microprocessor  software  are  presented  for  several 
computer  simulated  target  tracks. 
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I .  INTRODUCTION 


"In  December  1973  INTEL  Corporation  shipped  the  first 
8-bit  N-channel  microprocessor,  the  8080."  /~1_7  In  the 
months  that  followed,  several  semiconductor  manufacturers 
introduced  their  own  version  of  INTEL'S  8080  microprocessor 
The  resulting  competition  has  driven  microprocessor  prices 
down  and  performance  up.  Advanced  Micro  Devices  annoimced 
in  April  1976,  their  AM9080A-4  microprocessor,  with  an 
instruction  execution  time  one-half  that  of  INTEL'S 
original  8080,  was  available  for  $21.00  /~2_7.  In  the 
same  annoimcement ,  they  predicted  the  AM9080A  would  cost 
$5*00  by  1980.  Paralleling  the  advances  in  microprocessor 
technology,  the  manufacturers  of  semiconductor  memories  are 
making  steady  progress  in  expanding  the  memory  capability 
of  single  semiconductor  chips,  reducing  memory  access  time 
and  subsequently  reducing  the  cost  of  high-speed  semiconduc 
tor  memories  used  with  microprocessors. 

These  advances  in  semiconductor  technology  are  taking 
place  at  a  time  when  the  U.  S.  Navy  is  debating  the  cost 
of  rebuilding  the  fleet.  A  significant  portion  of  the  cost 
of  building  a  modern  warship  is  the  electronics  in  its 
weapons  systems.  Microprocessors  offer  the  potential  to 
(1)  reduce  the  cost  of  digital  systems,  (2)  perform  complex 
functions  at  remote  stations,  relieving  the  congestion  at 
larger  central  computing  facilities,  and  (3)  perform  func¬ 
tions  currently  handled  by  watch  personnel,  thus  reducing 
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the  manning  requirements  of  watch  sections. 


A.  OBJECTIVES  OF  THESIS 

The  primary  objective  of  this  thesis  is  to  demonstrate 
the  capability  of  microprocessors  to  perform  surface  search 
radar  tracking,  thus  relieving  the  tedium  of  this  function 
from  the  typical  underway  watch  section. 

The  secondary  objective  of  this  thesis  Is  to  investigate 
the  performance  of  a  Kalman  filter  in  the  two-dimensional 
polar  coordinate  system,  with  measurement  noise,  typical 
of  a  surface  search  radar. 

B.  SCOPE  OF  THE  WORK 

The  scope  of  the  work  was  limited  to  the  development  of 
a  Kalman  filter  tracking  algorithm  for  the  INTEL  8080 
microprocessor  and  simulation  of  circuitry  necessary  to 
interface  with  the  AN/SPS-10  surface  search  radar. 

A  Kalman  filter  was  chosen  because  it  is  computationally 
the  most  complex  of  the  tracking  filters  typically  used 
for  radar  tracking.  Thus,  it  would  be  the  most  demanding 
on  the  microprocessor.  Additionally,  the  Kalman  filter  is 
seldom  implemented  in  its  complete  form  due  to  its  computa¬ 
tional  complexity  thus  adding  additional  relevance  to  this 
investigation. 

The  parameters  for  the  AN/SPS-10  surface  search  radar 
were  used  for  computations  requiring  specific  radar  para¬ 
meters.  This  radar  was  selected  because  it  was  available 
at  the  Naval  Postgraduate  School  and  is  installed  aboard 
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the  majority  of  U.  S.  Naval  ships. 

INTEL  Corporation's  8080  microprocessor  was  used  for 
the  implementation  of  the  Kalman  filter  because  it  was 
available  at  the  Naval  Postgraduate  School  along  with  sup¬ 
porting  software. 

The  interface  between  the  microprocessor  and  the  radar 
was  not  constructed  due  to:  (1)  Cost,  time  and  manpower 
constraints,  and  (2)  The  performance  of  the  system  could 
not  be  effectively  measured  against  non-cooperative  targets 
of  opportunity. 

The  performance  of  the  radar  and  the  interface  equipment 
were  effectively  simulated  on  the  Naval  Postgraduate 
School's  IBM-360  computer,  thus  permitting  comparison  of 
the  microprocessor's  performance  against  that  of  the  IBM- 
360  system  for  known  target  tracks. 

The  output  of  the  tracking  system  would  be  displayed 
for  human  utilization  rather  than  being  eleotrically 
interfaced  to  fire  control  equipment.  The  minimum  standards 
of  performance  were  established  as  that  obtainable  from 
a  well  trained,  fully  manned  underway  watch  section  typical 
of  that  employed  on  a  modem  destroyer  not  equipped  with 
a  Naval  Tactical  Data  System  (NTDS). 
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II .  BACKGROUND 


A.  PERFORMANCE  OF  A  TYPICAL  CIC  WATCH  TEAM 

Onboard  naval  ships  not  equipped  with  Naval  Tactical 
Data  Systems  (NTDS),  the  tracking  of  radar  contacts,  both 
surface  and  airborne,  is  performed  manually  by  the  Combat 
Information  Center  (CIC)  watch  team.  During  normal  peace¬ 
time  steaming,  this  watch  team  may  consist  of  from  two  to 
ten  personnel,  depending  on  the  size  of  the  ship  and  the 
complexity  of  its  CIC  equipment.  Due  to  the  risk  of  col¬ 
lision  with  other  ships,  their  effort  is  primarily  devoted 
to  detecting,  tracking  and  computing  the  course,  speed, 
and  closest  point  of  approach  (CPA)  of  other  ships  within 
a  20  to  30  nautical  mile  radius.  This  can  be  a  tedious, 
error-prone  task  which  generally  requires  a  duplicate  effort 
by  the  bridge  watch  team  in  order  to  avoid  the  awesome 
consequences  of  a  collision  at  sea.  The  installation  of 
an  economical  device  to  track  and  compute  the  tactical 
data  for  surface  radar  contacts  could  eliminate  one  or  more 
watch-standers  in  CIC  while  improving  the  reliability  of 
the  information  provided  to  the  bridge  watch  team. 

The  manual  method  of  computing  the  course,  speed  and 
CPA  of  a  radar  target  involves  a  graphical  solution  of  the 
true  and  relative  velocity  vectors  of  the  radar  target's 
motion  and  own  ship's  course  and  speed.  The  target's 
relative  position  measurements  are  obtained  from  a  radar 


12 


repeater.  The  range  measurement  accuracy  is  at  "best  +  20 
yards  with  a  well  calibrated  repeater  and  a  well  trained 
operator.  The  bearing  measurement  accuracy  is  at  best 
+  1  degree  with  a  well  trained  operator.  The  own  ship's 
course  and  speed  used  in  determining  the  own  ship's  velocity 
vector  are  normally  based  on  the  ordered  course  and  speed. 
The  course  actually  steered  by  the  helmsman  may  vary  from 
the  ordered  course  by  +  2  degrees.  Depending  on  the  condi¬ 
tion  of  the  ship's  hull  and  its  engineering  plant,  the 
speed  will  be  accurate  within  +  1  knot.  The  combined 
results  of  these  errors  lead  to  a  target  course  accurate 
to  within  +  10  degrees  and  a  target  speed  accurate  to 
within  +  3  knots,  not  allowing  for  human  error.  The  major 
disadvantages  of  the  manual  solution,  in  addition  to  its 
poor  accuracy,  is  the  time  required  to  obtain  the  solution 
and  the  inability  to  obtain  solutions  if  either  the  target 
or  own  ship  maneuvers  between  radar  position  measurements, 
typically  3  or  more  minutes  apart.  The  estimates  of  error 
are  based  on  the  author's  personal  observation  of  unalerted 
CIC  watch  sections. 

B.  THE  TRACKING  PROBLEM 

The  "Tracking  Problem"  as  referred  to  in  this  investiga¬ 
tion  is  that  of  using  radar  measurements  to  accurately 
determine  the  course  and  speed  of  a  waterborne  radar  target. 
The  targets  of  interest  include  merchantmen,  warships,  and 
other  ships  capable  of  open  ocean  operations.  As  such. 
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their  motion  on  the  surface  of  the  ocean  can  he  approximated 
by  a  linear  dynamic  system.  In  a  sampled  data  system,  such 
as  a  search  radar  where  measurements  are  taken  once  each 
antenna  rotation,  the  target  dynamics  can  be  represented 
in  matrix  form  as  Z~3_7* 


Z 


n+1 


5  Z. 


n 


ra. 


n 


z 


n 


Y 


1  T  0  0 
0  10  0 
0  0  1  T 
0  0  0  1 


r 


tV2  0 

T  0 

0  T^/2 

0  T 


Z_  =  Vehicle  State  at  Scan  n 
n 

-  E/W  Postion 

» a 

-  E/W  Velocity 

-  N/S  Position 
Y^  -  N/S  Velocity 


T  =  Time  Between  Scan  n  and  n+1 


a  =  Acceleration  Acting  on  the  Target  from  Scan  n 
to  n  +  1  (E/W  and  N/S  components) 


If  the  radar  measurements  were  perfectly  accurate  and 
the  target  had  zero  acceleration,  then  the  tracking  problem 
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would  be  a  simple  velocity  vector  problem  (Figure  1). 
Unfortunately,  the  radar  bearing  and  range  to  the  target  on 
any  particular  scan  is  contaminated  by  errors  which  are  a 
function  of  the  radar  system  and  the  received  signal  to 
noise  ratio.  The  measurement  errors,  to  be  discussed  in 
the  next  section,  are  characteristically  Gaussian.  An 
optimal  estimator  is  well  suited  to  tracking  problems  such 
as  this . 

"An  optimal  estimator  is  a  computational  algorithm  that 
processes  measurements  to  deduce  a  minimum  error  estimate 
of  the  state  of  a  system  by  utilizing:  knowledge  of 
system  and  measurement  dynamics,  assumed  statistics  of 
system  noises  and  measurement  errors,  and  initial  condition 
information." 

A  Kalman  filter  is  a  common  optimal  filtering  technique 
for  estimating  the  state  of  a  linear  system  is 

chosen  for  implementation  because  it  is  computationally  the 
most  demanding  as  compared  to  other  techniques  used  in 
similar  tracking  systems  and,  as  such,  places  the  most 
demand  upon  the  microprocessor.  The  details  of  the  Kalman 
filter  will  be  covered  in  Chapter  IV  following  an  analysis 
of  the  measurements  and  associated  hardware  which  provide 
its  inputs  and  determine  its  parameters . 

C .  BtEASUREMENTS 

In  order  to  solve  the  tracking  problem,  it  is  necessary 
to  make  measurements  of  (1)  target  range,  (2)  target  bearing, 
(3)  own  ship's  course,  (4)  own  ship's  speed,  and  (5)  time 
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FIGURE  la  GEOGRAPHIC  MOTION  PLOT 
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Y 


FIGURE  Ic  VELOCITY  VECTOR  PLOT 
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of  the  measurement.  It  is  assumed  that  all  noise  is  Gaussian 
and  received  signal  to  noise  ratios  are  large  (»  1). 

The  target  range  is  determined  hy  measuring  the  time 
required  for  a  radar  pulse  to  travel  to  and  from  the 
target. 


R  -  Range  to  Target 
c  -  Speed  of  Light 

t  -  Elapsed  Time  from  Pulse  Transmission  to 
Pulse  Detection 

Using  threshold  detection  to  obtain  a  leading  edge  measure¬ 
ment,  the  theoretical  standard  deviation  of  error  for  the 
measurement  of  t  is  /~5_7* 

'  »  '[v 

Tr  (  ‘> 

- — I  ^  cW” 

(2  s/n)^  ^ 

At  =  Error  in  Elapsed  Time  Measurement  t 

=  Rise  Time  of  Received  Pulse 
r 

s/n  =  Signal  to  Noise  Ratio  of  Video  Pulse 

If  the  rise  time  of  the  pulse  (  Tj,)  is  limited  by  the  receiver 
IF  bandwidth  B,  then  ^  Z~5_7*  The  theoretical 

standard  deviation  of  range  error  using  leading  edge  detec¬ 
tion  is  therefore: 

_ c  ^ 

°R  ^  2B(2  s/n) 2 

cr^  =  Standard  Deviation  of  Range  Error 
The  mean  range  error  is  negligible  when  compared  to  the 
range  measurement.  Additionally,  the  mean  error  will  tend 


(At)^ 


18 


to  be  cancelled  by  subtraction  in  the  relative  velocity 
calculations .  Thus  the  computed  target  course  and  speed 
will  not  be  affected  by  the  mean  range  error. 

The  target  bearing  is  determined  by  measuring  the  bearing 
of  the  radar  beam  at  the  time  the  target  video  is  detected. 
The  mean  error  of  the  bearing  measurement  will  be  a  fimction 
of  the  antenna  beam  pattern,  and  the  target  range  and  radar 
cross  section.  The  resulting  bearing  to  target  will  not 
be  through  the  center  of  the  target  but  the  computed  target 
course  and  speed  will  not  be  affected. 

Using  beam  splitting,  a  technique  for  estimating  the 
center  bearing  of  the  target,  the  theoretical  minimum  stan¬ 
dard  deviation  of  bearing  error  for  a  Gaussian  beam  distri¬ 
bution  is  Z~5_7* 


where 


<^0 


1.06  0^ 

Ng(s/n)^ 


I®"''  U.C) 

V 


<f^  =  standard  Deviation  of  Bearing  Error 
0  =  Two  Way  Beamwidth  (=  0.4250  0-0=  Half 

g  D  D 

Power  Beamwidth)  (Same  units  as  Oq) 

(s/n)  =  Signal  to  Noise  Ratio  at  Center  of  Beam 

N  =  Number  of  pulses  transmitted  as  the  antenna 
°  rotates  through  2  0 

& 


The  standard  deviation  of  both  range  error  and  bearing 
error  are  functions  of  the  video  signal  to  noise  ratio. 

The  video  signal  to  noise  ratio  for  a  radar  with  single 
pulse  detection  is  given  by 
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n 


(47r)3  k  L^  R^ 

where 

Ppj^  =  Peak  Transmitted  Power  (watts) 

G  =  Antenna  Gain 
X  =  Wavelength  (meters) 

2 

O'  =  Radar  Cross  Section  of  Target  (m  ) 

ot  =  Attenuation  Constant  of  Propagation 
Medium  (db/m) 

R  =  Range  to  Target  (meters) 

k  =  Boltzmann's  Constant  (1.23  x  10~^^  joule/deg) 

T^  =  Standard  Temperature  (290°  K) 

=  Receiver  Noise  Bandwidth  (HZ) 

L  =  System  Losses 
s 

When  these  relationships  are  applied  to  the  AN/SPS-10 
(series)  Surface  Search  Radar,  with  the  AS-II6I/SPS 
Antenna,  the  standard  deviations  for  bearing  and  range  can 
be  related  to  target  range  and  cross  sectional  area  as 
shown  in  Figures  2  and  3*  The  minimum  performance  charac¬ 
ters  of  this  radar  system  are  listed  below. 


Peak  Power  =  I90  kw  Z”^_7 
Antenna  Gain  =  30  db  8_7 
Pulse  Repetition  Rate  =  625  HZ  Z~6_7 
Frequency  =  5825  MHZ  Z~^_7 
Propagation  Attenuation  =  .6  x  10“-^  db/m  /~5_7 
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FIGURE  2  STANDARD  DEVIATION  OF  RANGE  ERROR  VS.  RANGE 
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FIGURE  3  STANDARD  DEVIATION  OF  BEARING  ERROR  VS.  RANGE 


Noise  Figure  =  22.9  db 
Bandwidth  =  5  MHZ 

Noise  Bandwidth  (I.05  x  B)  Ci  5.25  MHZ 
System  Losses  5  db 
Antenna  Rotation  Rate  =  17  RPM 
Antenna  Horizontal  Beamwidth  =  1.9° 


^7_7 

Z-5J7 

/■5J7 

Zf8j7 


As  proposed  earlier,  a  threshold  detector  could  be  used 
to  indicate  the  occurrence  of  the  leading  edge  of  the 
received  radar  video  pulse.  The  threshold  to  noise  ratio 
of  such  a  detector  will  establish  the  ranges  at  which 
various  targets  can  be  tracked  and  the  minimxim  signal  to 
noise  ratio  which  the  tracking  system  will  encounter. 
Leading  edge  detection  is  preferred  because  it  eliminates 


errors  resulting  from  variations  in  the  pulse  width. 

When  considering  the  threshold  to  noise  ratio  of  such  a 
detection  scheme,  it  is  necessary  to  establish  the  desired 
mean  time  between  false  alarms,  and  the  desired  minimum 
probability  of  detection  of  the  target. 

The  mean  time  between  false  alarms  is  a  measure  of 
how  often  the  tracking  system  will  erroneously  recognize 
a  noise  pulse  as  a  target.  The  tracking  system  can  not 
function  effectively  if  this  occurs  frequently.  A  mean 
time  between  false  alarms  of  60  sec  has  been  selected.. 

The  mean  time  between  false  alarms  is  a  function  of  the 
signal  to  noise  ratio  and  the  IF  bandwidth  of  the  receiver 
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(-^Tp-^)  =  Threshold  to  Noise  Power  Ratio 


IF  “ 
FA  ^ 


IF  Bandwidth 

Mean  Time  Between  False  Alams 


The  exponential  tern  in  this  expression  is  the  inverse  of 
the  probability  of  false  alarm  Z~5_7« 


p  =  Exp(-  It  ) 

Tfo 


.  ,1 


c  .'b 


,i9 


Thus 


FA  T„^  B 


FA  SiF  -2 


The  mean  time  between  false  alams  as  seen  by  the  tracking 
system  can  be  increased  if  the  tracking  system  is  disabled 
whenever  targets  are  not  anticipated.  For  the  AN/SPS-10, 
the  maiximum  unambiguous  range  is  129.6  nm  of  which  only 
30  nm  is  of  interest  to  the  tracking  system.  Therefore,  ^ 
the  measurement  circuitry  could  be  disabled  77  percent  of 
the  time.  This  would  multiply  the  mean  time  between 
false  alarms  by  4.32.  An  additional  advantage  can  be 
gained  by  using  bearing  and  range  gates  to  further  limit 
the  opportunity  for  false  alarms.  If  such  bearing  and 
range  gates  isolated  1000  yard  squares  for  each  target 
and  the  tracking  system  were  tracking  eight  targets,  the 
mean  time  between  false  alarms  would  be  increased  by  an 
additional  factor  of  l4l4.  The  combination  of  these  factors 
leads  to  a  probability  of  false  alarm  of  2.036  x  10  ^  and 


i£. 


24 


a  threshold  to  noise  power  ratio  of  10.8  or  10. 3  db. 

The  probability  of  detection  of  the  target  is  a  measure 
of  how  often  the  tracking  system  will  recognize  the  target 
video  when  a  target  is  present.  It  is  essential  to  the 
effectiveness  of  the  tracking  system  that  it  detect  the 
target  video  reliably.  The  probability  of  detection  is  a 
complex  function  of  the  signal  to  noise  ratio,  the  charac¬ 
teristics  of  the  detector,  and  the  probability  of  false 
alarm  requiring  numerical  techniques  Z~5_7'  Figure  2,7 
of  Reference  5  shows  a  plot  of  the  probability  of  detection 
of  a  sine  wave  in  noise  as  a  function  of  the  signal  to 
noise  power  ratio  for  various  probabilities  of  false  alam. 
This  plot  shows  that  a  signal  to  noise  ratio  of  12.5  db  to 
14.0  db  is  required  for  a-  probability  of  detection  of  O.90 
to  0.99»  respectively,  and  a  probability  of  false  alarm 
of  2.0  X  10“-^ ,  If  all  targets  were  perfect  radar  reflec¬ 
tors,  it  would  be  possible  to  predict  the  performance  of 
the  threshold  detector  and  the  range  at  which  various 
targets  could  be  tracked.  However,  the  reflected  radar 
signals  from  the  targets  of  interest  will  fluctuate  due  to 
cancellation  and  reinforcement  of  the  waves  reflected 
from  various  portions  of  the  target.  These  fluctuations 
can  be  modeled  as  either  Swerling’s  Case  I  or  Swerling's 
Case  II  Z~5_7»  In  either  case,  an  additional  signal  to 
noise  ratio  of  8  db  to  17  db  will  be  required  to  insure 
a  probability  of  detection  of  0.9  to  0.99f  respectively 

zf5_7. 
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Thus  the  proposed  threshold  detection  scheme  should  have 


a  threshold  to  noise  power  ratio  of  10.8  to  achieve  an 
effective  mean  time  between  false  alarms  of  6o  sec.  This 
would  lead  to  a  worse  case  minimum  average  signal  to  noise 
ratio  of  20.5  db  to  31*0  db  for  a  probability  of  detection 
of  0.9  to  0.99 1  respectively. 

For  a  radar  measurement  system  with  the  parameters 
established  above,  targets  with  the  following  average 
radar  cross  sections  could  be  tracked  with  the  AN/SPS-10 
at  the  ranges  indicated. 


Average  Radar  Cross 
Section  (m^) 

1 

10 


10' 


10- 


10 


10- 


Range 

(yards/meters ) 

3,600/3,290 

6,400/5.850 

11,500/10,520 

20,000/18,290 

35,000/32,000 

60,000/54,860 


Based  on  a  threshold  to  noise  ratio  of  10.8,  the  stan¬ 
dard  deviation  of  the  radar  bearing  and  range  error  (Figures 
2  and  3)  will  be  O.0763  degrees  and  7*057  yards,  respec¬ 


tively. 

The  own  ship's  course  and  speed  measurements  come  from 
the  ship's  gyro  and  pitometer  log,  respectively.  Accordingly, 
their  error  statistics  are  a  function  of  the  equipment 
installed  on  the  particular  ship.  This  equipment  is  fre¬ 
quently  associated  with  the  ship's  fire  control  equipment 
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and,  as  such,  is  intended  for  tracking  operations. 

The  gyro  has  a  dual  influence  on  the  tracking  prohlem. 

It  not  only  determines  own  ship  course  hut  also  determines 
the  true  hearing  of  the  radar  aintenna.  The  mean  error  in 
radar  hearing,  and  likewise  own  ship's  course,  will  result 
in  an  error  in  the  solution.  This  can  he  viewed  as  a 
shift  in  the  reference  plane  which  does  not  prevent  the 
correct  prediction  of  the  future  relative  position  of 
targets.  As  a  refinement  to  the  tracking  system,  it  would 
he  better  to  measure  the  average  course  and  speed  made 
good  since  the  last  measurement  rather  than  the  instan¬ 
taneous  values  at  the  time  of  the  measurement.  Such  an 
average  measurement  could  he  achieved  hy  continuously 
integrating  the  variations  in  the  resolved  N/S  and  E/W 
velocities  using  a  time  constant  equal  to  one-half  the 
period  of  the  antenna  rotation  rate .  This  integration 
process  would  have  the  desirable  side-effect  of  minimizing 
the  variance  of  error  in  own  ship's  course  and  speed. 

For  subsequent  calculations,  it  is  assumed  that  the  standard 
deviation  of  error  of  own  ship's  gyro  is  0.1  degrees  and 
own  ship's  speed  indicator  is  0.1  knots. 

The  fifth  and  final  measurement,  the  time  of  the  measure¬ 
ment,  could  be  made  with  a  very  high  accuracy,  depending  on 
the  quality  of  the  clock  used.  The  objective  is  to  deter¬ 
mine  the  elapsed  time  between  position  measurements.  The 
antenna  rotation  rate  determines  the  mean  elapsed  time 
between  measurements  on  a  unique  target,  and  the  pulse 
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repetition  rate  of  the  radar  determines  the  useful  precision 
to  which  the  elapsed  time  needs  to  be  measured.  The  mean 
rotation  rate  for  the  AS-II61/SPS  Antenna  is  17  RPM.  The 
mean  pulse  repetition  rate  for  the  AN/SPS-10  (series) 

Radar  is  637 *5  HZ.  Thus,  the  mean  elapsed  time  between 
measurements  is  3*53  sec  and  the  useful  precision  of  the 
time  measurement  is  +  0.78^  ni  sec. 

D.  MICROPROCESSORS 

A  microprocessor  is  a  large  scale  integrated  (LSI) 
semiconductor  device  or  family  of  devices  which  contain 
an  arithmetic  logic  unit  (ALU)  and  several  working  registers 
capable  of  performing  binary  or  BCD  arithmetic  and/or 
Boolean  logic  operations.  Microprocessors  generally  fall 
into  one  of  five  categories;  (1)  single-chip  controllers, 

(2)  4- bit  family,  (3)  8-bit  family,  (4)  l6-bit  family, 
and  (5)  bit  slices  /~9_7* 

The  single  chip  controllers  are  LSI  devices  which 
incorporate  all  the  necessary  hardware,  including  program 
and  data  memory  circuitry,  on  a  single  chip.  Their  data 
processing  applications  include  devices  such  as  calculators, 
gas  pumps,  cash  registers,  and  scales.  They  include  devices 
with  4-bit  and  8-bit  ALUs.  Typically  the  4-bit  versions 
perform  serial  BCD  operations.  They  are  limited  in  their 
speed  of  program  execution  by  the  power  of  their  instruc¬ 
tion  sets  and  the  size  of  their  ALUs.  The  complexity  of 
the  operations  they  can  perform  is  limited  by  the  size  of 
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their  memories.  Some  manufacturers  are  developing  single¬ 
chip  versions  of  their  8-bit  family  devices  which  will 
expand  the  capabilities  of  the  single-chip  controllers 
to  the  point  of  overlapping  the  8-bit  families 

The  if— bit  microprocessor  families  are  similar  in 
function  to  the  single-chip  controllers  except  the  memory  is 
not  included  on  the  chip.  Therefore,  the  4-bit  families 
are  capable  of  handling  larger  programs.  The  term  family 
refers  to  the  set  of  compatible  chips  which  go  together 
to  make  up  the  microprocessor.  They  require  greater  design 
effort  than  the  single  chip  controllers  because  several 
chips,  including  clocks  and  l/O  interface  chips,  along 
with  the  separate  memory,  are  required  for  even  the  simplest 
task.  Like  the  4-bit  controllers,  the  4-bit  microprocessor 
families  are  limited  in  the  speed  of  their  program  execution 
by  the  power  of  their  instruction  sets  and  the  size  of 
their  ALUs. 

The  8-bit  families  are  presently  dominating  the  micro¬ 
processor  applications  Z~9_7*  Their  instruction  sets  are 
approaching  the  complexity  of  those  used  in  minicomputers. 
They  typically  provide  for  8-bit  data  manipulation  along 

with  l6-bit  memory  addressing.  " - INTEL  Corporation's 

8080  family,  with  its  enhanced  8O8OA  CPU,  Motorola  Semi¬ 
conductor's  6800  family,  with  its  enhanced  6800D  CPU, 
and  Rockwell's  PPS-8  family  currently  rank  one,  two,  and 
three  in  popularity  among  users."  Z~9_7 

The  strength  of  the  8O8O  lies  in  the  power  of  its 


29 


instruction  set  and  the  extent  of  its  software  support. 

The  8080  has  72  basic  instructions,  including  conditional 
branching,  decimal  as  well  as  binary  arithmetic,  logical, 
register  to  register,  stack  control  and  memory  reference 
instructions  These  instructions  can  be  executed 

in  4  to  18  cycles  of  its  2  MHZ  clock  allowing  some  instruc¬ 
tions  to  be  executed  in  2  psec  1_7«  These  instructions 
facilitate  l6-bit  addition  and  register  to  register  trans¬ 
fers.  Among  its  software  support  are  a  System  Monitor,  an 
Assembler,  a  Text  Editor,  and  a  PL/M  compiler,  cross 
assembler,  and  simulator.  The  popularity  of  PL/M,  an 
INTEL  derivative  of  PL/l,  and  the  8080  have  caused  other 
manufacturers  to  model  their  microprocessors  after  the 
8080,  such  as  Advanced  Micro-Devices  9080A,  or  to  develop 
their  own  version  of  PL/M  such  as  PL/M  6800  by  Intermetrics 
for  the  Motorola's  6800. 

The  8-bit  devices  are  limited  in  performance  by  their 
8-bit  data  path  to  and  from  memory.  This  necessitates 
using  several  memory  cycles  to  read  or  write  address  and 
other  l6-bit  variables. 

The  l6-bit  microprocessor  family,  including  National's 
PACE,  Texas  Instruments  Inc. 's  9900,  and  General  Instrument 
Corp. 's  CPI600,  attempt  to  reduce  the  number  of  instructions 
and  machine  cycles  required  to  manipulate  memory  addresses 
(typically  l6-bit  variables),  and  data  originating  from 
A-D  converters  or  going  to  D-A  converters  (typically  10-16 
bits).  Their  performance  approaches  that  of  minicomputers. 
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Western  Corp.'s  MCP  l600  was  originally  designed  for  DEC'S 
LSI-11  l6-bit  minicomputer  and  was  subsequently  distributed 
as  an  independent  device  /~9_7* 

The  bit-slice  architectures,  such  as  Advanced  Micro 

Devices  Inc.'s  AM  2900,  INTEL  Corp.'s  3002,  and  Texas 

Instrument  Inc.'s  SBP0400,  are  taking  on  the  tough  process- 

control  and  high  speed  controller  jobs  previously  handled 

by  minicomputers.  A  typical  bit-slice  device  is  a  2  or  4 

bit  vertical  slice  of  a  CPU  which  is  horizontally  expandable 

in  2  or  4  bit  increments  to  form  a  16  bit,  32  bit  or  other 

multiple  of  2  or  4  bit  CPU.  They  utilize  Schottky  TTL, 

2 

I  L  or  ECL  circuitry  and  thus  are  much  faster  than  the 
previously  mentioned  microprocessors  which  are  primarily 
MOS  devices.  The  bit-slice  architecture  offers  the  ulti¬ 
mate  in  design  flexibility  while  demanding  the  most  effort 
by  the  designer.  The  design  effort  may  include  an  extensive 
amount  of  microprogramming  prior  to  development  of  the 
working  program.  Additionally,  bit-slice  devices  are  the 
most  expensive  form  of  the  microprocessor. 

In  summary,  the  8-bit  and  16-bit  microprocessors  are 
the  best  alternatives  to  the  tracking  problem.  The  single¬ 
chip  controllers  and  4-bit  families  have  insufficient  memory 
and/or  are  too  slow  to  handle  the  Kalman  filter  algorithm. 
The  bit  slice  architectures  have  excess  capability  while 
requiring  unnecessary  design  and  hardware  cost.  As  will 
be  seen  later,  the  tracking  problem  can  be  solved  accurately 
with  l6-bit  resolution  in  the  arithmetic  operations.  Thus 


31 


an  8-bit  microprocessor  programmed  for  double  precision 
arithmetic  or  a  l6-bit  microprocessor  program  for  single 
precision  arithmetic  is  suitable. 
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Ill .  SYSTEM  INTERFACES 


The  required  interface  equipment  can  he  grouped  into 
three  categories;  (1)  the  radar  interface,  (2)  the  own 
ship  course/speed  interface,  sind  (3)  the  operator/user 
interfaces.  Although  the  interface  hardware  was  not  con¬ 
structed,  it  is  essential  to  the  validity  of  the  proposed 
tracking  system  to  verify  the  feasibility  of  making  the 
required  radar  and  own  ship  measurements  within  the  pro¬ 
posed  tolerances. 

A.  RADAR  INTERFACE 

The  radar  interface  serves  four  functions;  (1)  to  con¬ 
vert  the  radar  video  to  rectangular  pulses  compatible 
with  the  logic  circuitry  in  the  range  measurement  subsystem, 
(2)  to  identify  the  radar  video  of  the  target  whose  position 
is  to  be  measured,  (3")  to  measure  and  store  the  bearing, 
range  and  time  of  measurement  of  the  designated  radar 
target,  and  to  transmit  a  signal  to  the  own  ship 

course/speed  and  the  operator/user  interfaces  to  indicate 
when  the  radar  measurements  have  occurred. 

The  conversion  of  the  radar  video  to  a  rectangular  pulse 
compatible  with  the  bearing  and  range  measurement  circuitry 
could  best  be  accomplished  at  the  threshold  detector.  The 
threshold  detector  could  consist  of  one  of  a  variety  of  LSI 
voltage  comparators  currently  available .  The  input  video 
would  be  in  the  0-2  volt  range  depending  on  the  setting  of 
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the  radar  IF  aunplifier  gain  adjustment.  The  threshold 
voltage  would  he  a  technician  adjustment  made  hy  first 
setting  the  IF  amplifier  gain  to  obtain  a  good  video 
display  on  the  radar's  master  repeater  and  then  adjusting 
the  threshold  voltage  on  the  comparator  to  just  detect  a 
known  signal.  The  critical  parameters  of  the  threshold 
detector  are  a  bandwidth  greater  than  that  of  the  IF 
amplifier  (5  MHZ)  and  a  high  noise  immimity. 

Target  identification  could  be  accomplished  either 
with  software  or  hardware .  A  software  solution  would  result 
in  the  measurement  and  processing  of  a  large  number  of 
targets  which  the  operator  had  no  interest  in  and  would 
place  an  unnecessary  burden  on  the  microprocessor.  Thus 
a  hardware  technique  is  preferred.  A  simple  solution  to 
this  problem  is  the  use  of  bearing  and  range  gates  to 
enable  the  measurement  circuitry  at  the  time  a  particular 
target  video  is  anticipated. 

The  width  and  time  of  occurrence  of  the  bearing  and 
range  gates  can  be  determined  in  several  ways,  all  of  which 
are  a  function  of  the  accuracy  with  which  the  future  loca¬ 
tion  of  the  particular  target  can  be  estimated.  This 
estimate  is  a  function  of  (1)  the  maneuverability  of  the 
target,  (2)  the  elapsed  time  between  measurements  and  (3) 
the  accuracy  of  the  computed  course  and  speed  of  the  target. 

The  majority  of  targets  of  interest  will  be  capable  of 

turning  radii  of  100-1000  yards,  velocities  of  0-40  knots, 

2 

and  accelerations  of  0-.5  yards/sec  .  The  antenna  rotation 
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rate  determines  the  mean  time  between  measurements ,  which 
in  the  case  of  the  AS-II6I/SPS  antenna  equates  to  3*53 
seconds  Z~8_7.  The  accuracy  of  the  computed  course  and 
speed  of  the  target  will  depend  on  the  tracking  system 
parameters  and  the  maneuverability  of  the  target. 

After  the  three  major  factors  have  been  accoimted  for, 
the  probable  measurement  error  should  be  added  .to  accoxmt 
for  differences  in  the  target’s  predicted  bearing  and  range 
from  the  measurement  system  bearing  and  range  values. 
Because  of  the  uncertainty  of  these  many  factors,  it  is 
necessary  to  use  wide  bearing  and  range  gates  to  insure 
that  the  target  is  not  lost.  The  \mfortunate  outcome  of 
using  wide  bearing  and  range  gates  is  the  increased  likeli¬ 
hood  that  two  adjacent  targets  will  enter  each  others 
bearing  and  range  gates  and  become  confused. 

The  best  solution  to  this  dilemma  is  to  allow  the 
operator  to  select  a  target  classification  from  his  control 
station  which  would  determine  which  of  several  predeter¬ 
mined  bearing  and  range  gate  widths  or  width  determining 
algorithms  would  be  used  with  the  designated  target.  As 
will  be  shown  later,  the  designation  of  target  maneuver¬ 
ability  by  the  operator  is  not  only  useful  in  the  measure¬ 
ment  circuitry  but  also  enhances  the  Kalman  filter. 

The  determining  factor  in  the  accuracy  to  which  the 
bearing,  range  and  time  measurements  can  be  represented 
in  the  microprocessor  is  the  value  represented  by  the 
least  significant  bit  in  the  binary  data  word.  If  the 
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digitizing  process  is  carried  out  with  the  digital  zero 
equal  to  the  analog  zero  it  will  introduce  a  mean  error 
equal  to  one-half  of  the  value  represented  by  the  least- 
significant  bit  and  a  standard  deviation  of  error  equal  to 
3”^  times  the  value  represented  by  the  least  significant 
bit.  If  the  digital  zero  is  offset  by  minus  the  previous 
mean  error,  then  the  mean  error  will  be  zero  and  the 

_  i 

standard  deviation  of  error  equal  to  12~®  times  the  value 
represented  by  the  least  significant  bit.  Obviously,  the 
zero  offset  approach  is  preferred. 

The  length  of  the  binary  data  word  used  to  represent 
the  bearing,  range,  and  time  measurements  is  determined 
by  the  value  represented  by  the  least  significant  bit, 
and  the  maximum  value  which  is  to  be  represented. 

For  the  system  under  consideration,  the  least  signifi¬ 
cant  bit  should  represent  the  standard  deviation  of  error 
of  the  target  bearing  with  gyro  error  (O.I258  d’egrees  or 
0.002196  radians)  and  the  target  range  (7.057^  yards), 
and  the  useful  precision  of  the  time  measurement  (+  0.784 
msec),  as  a  minimum.  The  maximum  value  of  each  measurement 
should  be  2 ‘Tf  radians,  60,000  yards  and  approximately  8.0 
secs.  This  would  require  11.5  bits  for  the  bearing  measure¬ 
ment,  13.1  bits  for  the  range  measurement  and  13*3  bits 
for  the  time  measurement.  If  an  8-bit  microprocessor  is 
to  be  used,  then  data  words  of  8,  16,  or  24  bit  lengths 
are  practical.  In  addition,  the  INTEL  8O8O  microprocessor 
has  some  instructions  for  manipulating  l6-bit  data  words. 
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Thus  16-bit  data  words  were  considered  optimum.  This 
results  in  the  following  data  format: 


Measurement 

Bearing 

Range 

Time 


Range 
0-2-tf  rad 
0-60,000  yds 
0-8.0  sec 


MSB 
Tf  rad 

30,000  yds 
4.0  sec 


LSB 

0.0959  mrad 

0.9155  yds 

0.1221  msec 


By  changing  the  most  significant  bit  to  the  next  even  power 
of  two,  the  units  can  be  rationalized  as  follows: 


Measurement 

Bearing 

Range 

Time 


Range 
O-Btt  rad 

0-65,535  yds 

0-7«999  sec 


MSB 

4 . 0  rad 
32,768  yds 
4.0  sec 


LSB 

2"^^  rad 


1 . 0  yds 
2”^^  sec 


The  resolution  resulting  from  this  format  will  result  in 
less  than  1.0  percent  increase  in  the  standard  deviation 
of  error  as  a  result  of  digitizing  the  data.  The  remaining 
question  is  what  resolution  is  possible  with  the  available 
hardware . 

The  bearing  measurement  could  be  conveniently  made 
by  combining  one  of  the  synchro-to-digital  converter  modules 
such  as  Transmagnetics  series  I65I,  which  has  l4-bit 
resolution,  and  two  l6-bit  latchs ,  enabled  by  the  video 
circuitry,  to  store  the  first  and  last  bearings  from 
which  video  was  detected  from  a  particular  target  on 
each  scan.  The  microprocessor  software  could  then  compute 
the  center  bearing  to  the  target.  An  alternate  scheme 
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would  "be  to  use  a  single  latch  to  store  the  first  hearing 
from  which  video  was  detected  from  a  particular  target  on 
each  scan  and  a  counter  to  count  the  number  of  video  pulses 
from  the  target  on  each  scan.  The  microprocessor  software 
could  then  use  this  data  to  estimate  the  center  bearing. 

A  third  and  far  less  accurate  scheme  would  be  to  simply 
use  the  bearing  of  the  first  video  pulse.  The  first 
scheme  is  recommended. 

The  rainge  measurement  is  essentially  a  time  interval 
measurement.  A  range  of  1  yard  corresponds  to  6.1  nsec. 
Thus  a  binary  counter  would  have  to  be  clocked  at  163*9 
MHZ  to  count  radar  range  in  intervals  of  1  yard.  Ordinary 
Transistor-Transistor-Logic  (TTL)  counters  are  limited  to 
clock  frequencies  of  30  MHZ  and  below,  while  high-speed  TTL 
counters  can  operate  as  high  as  50  MHZ  Z~10_7.  Schottky- 
Clamped  TTL  counters  can  operate  reliably  at  clock  fre¬ 
quencies  of  100  MHZ  Z~10_7.  Emitter-Coupled-Logic  (ECL) 
circuitry  could  perform  at  higher  clock  frequencies  but 
the  complexity  and  cost  of  the  circuitry  would  increase 
considerably  due  to  the  unavailability  of  ECL  Medium-Scale 
Integrated  (MSI)  semiconductor  chips  with  the  desired 
functions.  A  Shottky-C lamped  TTL  Counter,  such  as  the 
74S197»  is  suggested,  with  a  compromise  reduction  in  the 
clock  frequency  to  81. 96  IVIHZ  resulting  in  a  minimum  range 
resolution  of  2  yards.  An  alternate  scheme  would  be  to 
mix  ECL  circuitry  with  the  Shottky  counters. 

Another  consideration  in  the  range  measurement  counter 
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is  the  propagation  of  carry  pulses  through  the  counter.  The 
7^3197  binary  counter  has  a  maximum  propagation  time  of  13 
nsec  through  each  stage.  For  the  15  stages  required  for  a 
maximum  range  of  65 » 536  yards ,  a  total  propagation  time  of 
195  nsec  would  be  required  for  a  carry  propagation  to  the 
highest  significant  bit.  This  is  16  times  the  clock  period, 
which  would  lead  to  unstable  data  at  the  counter  outputs. 

The  easiest  solution  to  this  dilemma  is  to  stop  the  clock 
input  to  the  counter  at  the  time  the  target  video  is 
detected  and  store  the  counter  value  after  it  has  stabilized. 
This  would  necessitate  a  separate  range  coimter  in  the 
interface  hardware  for  each  target  being  tracked. 

The  time  measurement  interface  is  comparatively  straight 
forward.  A  l6-bit  counter  driven  by  a  8,192  HZ  clock,  with 
a  storage  latch  for  each  target  will  suffice.  As  in  the 
range  counter,  the  carry  propagation  delay  must  be  taken 
into  account  when  storing  the  time  measurement.  Due  to 
the  slower  clock  frequency,  the  counter  output  will  be 
stable  99  percent  of  the  time.  Simple  combinational  logic 
circuitry  could  be  employed  to  synchronize  the  time  measure¬ 
ment  with  the  clock  to  avoid  storing  unstable  counter 
outputs . 

Using  currently  available  hardware,  the  radar  interface 
can  be  accomplished  with  resolutions  well  within  the 
radar  measurement  errors  specified  earlier  in  this  section. 
The  following  data  formats  are  recommended. 
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Measurement 

Ran^ 

No.  of  Bits 

Resolution 

Bearing 

Q-2-ir  rad 

14 

•3835  mrad 

Range 

0-65,535  yards 

15 

2 . 0  yards 

Time 

0-8  sec 

16 

.1221  msec 

B .  OWN  SHIP 

COURSE  AND  SPEED 

INTERFACE 

The  own  ship  course  and  speed  interface  serves  two 
fimctions;  (1)  to  smooth  the  variations  in  own  ship's 
course  and  speed,  and  (2)  to  measure  and  store  the  smoothed 
course  and  speed  to  be  associated  with  a  particular  target. 
There  are  two  options  of  how  to  perform  these  functions, 
stemming  from  the  fact  that  the  microprocessor  software 
uses  own  ship's  course  and  speed  data  solely  to  compute 
own  ship  north/south  and  east/west  velocities. 

The  first  approach  would  be  to  process  the  course  and 
speed  measurements  independently.  The  course  measurement 
could  be  smoothed  either  mechanically  or  with  an  over-damped 
servo  feedback  system.  A  synchro  to  digital  angle  converter 
could  be  employed  to  digitize  the  measurement.  The  speed 
measurement  could  be  smoothed  electrically  and  digitized 
using  an  analog  to  digital  converter. 

The  second  approach  would  be  to  combine  the  course  and 
speed  measurements  into  north/south  and  east/west  velo¬ 
cities.  The  course  synchro  signal  could  be  converted  to 
two  analog  signals  corresponding  to  the  sine  and  cosine  of 
own  ship's  course  using  a  synchro  to  DC  sine/cosine  converter, 
such  as  Transmagnetics  model  655N  V  (series).  These  analog 
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signals  could  then  he  multiplied  hy  an  analog  speed  signal 
with  two  DC  analog  multipliers.  Low-pass  filters  could  be 
used  to  smooth  the  two  signals  and  then  digital  to  analog 
converters  used  to  digitize  them.  This  approach  would 
simplify  the  work  of  the  microprocessor  but  additional 
errors  would  be  generated  in  the  analog  multipliers . 

A  compromise  approach  would  be  to  use  the  first  system 
but  substitue  a  synchro  to  digital  sine/cosine  converter 
in  place  of  the  digital  angle  converter.  With  this  tech¬ 
nique,  the  microprocessor  would  receive  smoothed  north/ 
south  and  east/west  course  component  measurements  and 
smoothed  speed  measurements . 

Assuming  a  standard  deviation  of  error  in  the  course 
measurement  of  0.1  degree  (1.745  x  10~^  radians)  prior 
to  digitizing,  a  digital  data  format  of  11.8  bits  would  be 
required  for  minimum  resolution.  Using  a  data  format 
similar  to  that  used  for  the  target  bearing  (MSB  =  2.0 
radians)  and  l4-bit  resolution,  the  standard  deviation 
of  error  will  only  increase  0.3  percent  as  a  result  of 
digitizing.  Assuming  a  standard  deviation  of  error  in  the 
speed  measurement  of  0.1  knot  prior  to  digitizing,  and  a 
maximum  speed  of  100  knots,  a  digital  data  format  of  10 
bits  would  be  required.  Using  a  12-bit  analog  to  digital 
converter  or  a  12-bit  synchro  to  digital  converter  with  the 
most  significant  bit  representing  64  knots,  the  standard 
deviation  of  error  will  increase  only  1.6  percent,  as  a 
result  of  digitizing. 
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In  summary,  the  own  ship  course  and  speed  interface 
could  he  accomplished  using  available  hardware,  with  the 
following  data  fomats  and  measurement  precision. 


Measurement 

Range 

No.  of  Bits 

Resolution 

Course 

0-2'rtrad 

14 

.3835  mrad 

Speed 

0-128  knots 

12 

.03125  knots 

C .  OPERATOR/USER  INTERFACE 

The  operator/user  interface  serves  two  functions;  (1) 
to  provide  operator  control  of  the  tracking  system,  sind 
(2)  to  display  the  output  data  to  the  users. 


The  operator  should  have  control  over  (1)  which  radar 
video  is  to  be  tracked,  (2)  the  number  of  targets  being 
tracked,  and  (3)  the  tracking  parameters  to  be  used  with 
each  target. 

It  is  proposed  that  the  operator  would  control  the 
system  in  the  following  manner.  The  measurement  interface 
circuitry  would  be  divided  into  tracking  elements,  each 
capable  of  handling  one  target.  Each  element  would  be 
composed  of  a  bearing  and  range  gate  generation  circuit, 
a  range  measurement  circuit,  and  storage  circuits  for  the 
bearing,  range,  time,  own  ship's  course  and  own  ship's 
speed  data. 

The  operator's  control  panel  would  have  an  off/on/start 
switch  for  each  tracking  element.  The  system  would  be 
turned  on  with  all  tracking  elements  turned  off .  To  track 
a  target,  the  operator  would  place  his  radar  repeater's 
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bearing  and  range  cursor  on  the  target  video  to  be  tracked 
and  momentarily  place  the  tracking  element  switch  in  the 
start  position  when  the  radar  sweep  passed  the  target  and 
then  place  the  tracking  element  switch  in  the  on  position. 
This  process  would  momentarily  couple  the  radar  repeater's 
bearing  and  range  cursor  video  to  the  tracking  elements 
bearing  and  range  gate  circuit  and  enable  the  tracker's 
measurement  circuitry,  thus  taking  a  measurement  of  the 
target's  position.  If  the  measurement  circuitry  detected 
target  video,  an  indicator  light  would  be  lit.  To  stop 
tracking  a  target,  the  operator  would  simply  turn  the 
tracking  element  off.  A  control  data  word  (8  bits)  would 
be  periodically  read  by  the  microprocessor  to  advise  it  as 
to  which  tracking  elements  have  updated  measurement  data 
available.  The  operator  would  have  a  second  multi-position 
switch  which  would  allow  him  to  select  the  bearing  and 
range  gate  widths  and  the  variance  of  acceleration  para¬ 
meters  to  be  used  with  each  target. 

The  output  data  to  be  displayed  would  include  (1)  target 
designation,  (2)  target  bearing,  (3)  target  range,  (4)  tar¬ 
get  course,  and  (5)  target  speed,  as  a  minimum.  For  more  . 
sophisticated  systems  the  target's  closest  point  of  approach 
(CPA)  could  be  computed  and  displayed  as  (1)  CPA  bearing, 

(2)  CPA  range  and  (3)  time  of  CPA.  The  speed  of  the  data 
output  is  critical  to  the  capability  of  the  system  to  track 
multiple  targets.  Direct  memory  access  techniques  are 
suggested  as  the  only  practical  solution. 
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After  the  display  data  is  extracted  from  the  tracking 
system,  several  display  techniques  could  be  used.  Tabular 
displays  in  the  form  of  7-segment  LED  modules  or  cathode 
ray  tubes  would  facilitate  displaying  the  data  to  several 
users  at  the  operators  station  and  also  transmission  of 
the  data  to  remote  stations,  such  as  the  bridge,  CO's 
cabin,  and  weapons  stations.  One  of  the  major  advantages 
of  an  automatic  tracking  system  is  its  capability  to 
instantly  and  continuously  disseminate  tactical  data  to 
remote  stations  thus  eliminating  the  need  for  passing 
this  information  on  sound  powered  phones  as  is  currently 
the  practice  on  non-NTDS  ships . 
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IV.  KALMAN  FILTER 


As  stated  earlier,  the  Kalman  filter  is  a  common  optimal 
filtering  technique  for  estimating  the  state  of  a  linear 
system.  "Among  the  presumed  advantages  of  this  type  of 
data  processor  are  that  it  minimizes  the  estimation  error 
in  a  well  defined  statistical  sense  and  that  it  utilizes 
all  measurement  data  plus  prior  knowledge  about  the  system. 
The  corresponding  potential  disadvantages  are  its  sensitiv¬ 
ity  to  erroneous  a  priori  models  and  statistics,  and  the 
inherent  computational  burden."  _/ 

A.  DISCRETE  KALMAN  FILTER  ALGORITHM 

The  rotation  of  the  surface  search  radar  antenna 
results  in  periodic  measurements  of  the  target  positions. 
Accordingly,  the  discrete  form  of  the  Kalman  filter,  rather 
than  the  continuous  form,  is  appropriate.  The  discrete 
Kalman  filter  can  be  described  by  a  set  of  matrix 
equations . 

Using  the  target  state  (Z^)  representation  presented  in 
the  description  of  the  tracking  problem,  in  Section  I  B, 
the  matrix  form  of  the  measurement  process  is  /~3_75 


w  =  HZ  +  V 
n  n  n 


w  = 
n 


‘*5 


Xjjj(n)  =  Measured  x  Coordinate  at  Scan  n 
y-j^Cn)  =  Measured  y  Coordinate  at  Scam  n 


H  = 


10  0  0 
0  0  10 


=  Target  State  Vector  at  Scan  n 


V  = 

n 


v^(n) 

Vy(n) 


V  (n)  =  Random  Measurement  Error  on  x  (n) 

X.  in 

Vy(n)  =  Random  Measurement  Error  on 
The  estimated  state  vector  at  scan  n  (Z  )  is 

\  -  S  ^n-1 

i  -  (See  Tracking  Problem  Page  14) 

Z  ^  =  Optimal  Estimate  of  Target  State  Vector 
at  Scan  n-1 


The  optimal  estimate  of  the  target  state  at  scan  n  (Z^) 


is  given  by  Z~3_7! 


n 


Kn  = 


,  +  K 

n  n 

( w 
n 

■  H 

’xx 

V 

!  /T 

XX  ' 

®xy 

/T 

‘yx 
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!  /T 

yx  ' 

®yy 

/T 
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=  Gain  Matrix  at  Scan  n 


The  gain  matrix  (K^)  is  computed  from  /~3_7* 


rn  rp  -t 

K  =  P  H^(HP  +  R  )  ^ 
n  n  '  n  n^ 
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P  =  Estimated  Covariance  of  Estimation  Errors 
^  Matrix  Prior  to  Processing  Scan  n  Measurements 


a/(n) 


o-/(n) 


=  Covariance  of  Measurement 
Error  Matrix  at  Scan  n 


Following  the  processing  of  scan  n  measurements,  an 
optimal  estimate  of  the  covariance  of  estimation  error 
matrix  for  scan  n  can  be  computed  Z”3_7** 

P  =  (I  -  K  H)P 
n  n  n 

I  =  Identity  Matrix 


From  this,  it  is  possible  to  compute  the  estimate  of 
the  covariance  of  estimation  error  matrix  for  scan  n  +  1 

rsy. 

^n+l  "  *n+l  *n+l  ^  n+1  ®  „+! 

2 

Q  ~  cr  =  Variance  of  Random  Target  Acceleration 

3. 


The  basis  for  the  discrete  Kalman  filter  in  this  form 
is  the  legitimacy  of  the  target  maneuver  model.  The 
variance  of  target  acceleration  must  be  equal  and  indepen¬ 
dent  along  the  x  and  y  axis,  and  from  scan  to  scan.  This 
is  true  for  the  targets  of  interest.  Thus,  this  fora  of 
the  discrete  Kalman  filter  is  valid  if  the  position 
measurements  were  made  in  the  x  and  y  coordinates. 

Using  the  polar  coordinate  system  shown  in  Figure  4, 
the  following  transformations  can  be  used  with  the  radar's 
bearing  (9)  and  range  measurements  (R). 
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Y 


FIGURE  4  COORDINATE  SYSTEM 


48 


X  =  R  sin  0 
m 

ym  =  R  cos  0 

Using  these  transformations,  the  elements  of  the  covariance 
of  measurement  error  matrix  are  Z~ll_7: 


<7-2 

r 

sin^0 

+  R^ 

°0 

cos  0 

cos^0 

+  R^ 

^  2 
^0 

.  2 
sin 

i  sin  2  0 

2 

cr^  =  Variance  of  Range  Error 
2 

o'q  =  Variance  of  Bearing  Error 

This  modification  completes  the  description  of  the  discrete 
Kalman  filter. 

In  some  applications,  it  is  possible  to  precompute  the 
values  of  the  gain  matrix  and  store  them  in  the  computers 
memory  for  subsequent  use.  This  greatly  simplifies  the 
software  algorithm.  There  are  two  reasons  why  this  is  not 
possible  in  this  case. 

The  first  is  that  the  elapsed  time  between  measurements 
(T)  is  not  constant.  This  is  primarily  due  to  the  fact  that 
the  radar  antenna  rotates  relative  to  own  ship’s  heading. 
Thus,  changes  in  own  ship's  course  can  significantly 
increase  or  decrease  the  elapsed  time  between  measurements 
(T). 

The  second  reason  is  that  the  covariances  of  measurement 
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error  (cr^,  (Ty,  o^y)  are  a  function  of  the  hearing  and  range 
to  the  target.  Consequently,  the  gain  matrix  must  be  com¬ 
puted  on  line  for  each  target  being  tracked.  A  flow  chart 
of  these  calculations  is  shown  in  Figure  5* 

The  solution  of  the  Kalman  filter  equations  in  matrix 
form  is  unnecessarily  complex  for  the  microprocessor. 
Therefore,  an  algebraic  solution  for  the  value  of  each 
element  in  the  various  matrices  was  formulated.  Appendix 
A  lists  this  algebraic  form  of  the  Kalman  filter  equations. 

B.  KALMAN  FILTER  PARAMETERS 

The  performance  of  the  Kalman  filter  is  determined  by 

three  parameters?  (1)  variance  of  bearing  error  (Oq  ), 

2 

(2)  variance  of  range  error  (cr'^  ),  and  (3)  variance  of 

2 

target  acceleration  {cr  ).  The  first  two  parameters 
are  functions  of  the  radar,  the  ship's  gyro,  and  the  inter¬ 
face  system.’  The  third  parameter  is  a  function  of  the 
target's  relative  motion  and  the  desired  smoothness  of  the 
output  data. 

The  variance  of  bearing  error  is  equal  to  the  sum  of 
the  variance  of  the  radar  bearing  error,  (.0058  deg  ),  the 
variance  of  the  bearing- interface  digitizing  error  (.00016 
deg  ) ,  and  the  variance  of  the  ship's  gyro  error  (0.01  deg  ) . 
Thus,  the  total  variance  of  bearing  error  is  0.016  deg  or 
.0000049  rad^. 

The  variance  of  range  error  is  equal  to  the  si^m  of  the 

2  2 

variance  of  the  radar  range  error  (50  yd  or  42  m  )  and 


50 


51 


FIGURE  5  DISCRETE  KALMAN  FILTER  FLOW  DIAGRAM 


the  variance  of  the  range  interface  digitizing  error 
2  2 

(1.3  yd  or  1.1  m).  Thus  the  total  variance  of  range 

2  2 

error  is  equal  to  5i  yd  (43  m  ) . 

The  variance  of  target  acceleration  parameter  is  a 
compromise  between  the  user's  desired  smoothness  of  the 
output  data  for  non-maneuvering  targets  and  the  ability  of 
the  filter  to  provide  the  closest  estimate  of  a  target's 
position  and  velocity  for  a  maneuvering  target. 

A  possible  solution  to  this  dilemma  is  to  allow  the 
operator  to  select  the  variance  of  acceleration  parameter 
to  be  used  with  each  target.  This  parameter,  along  with 
the  bearing  and  range  gate  widths,  could  be  selected  by 
enabling  one  of  several  read-only-memories  (ROMs)  containing 
the  desired  parameters.  The  particular  ROMs  installed 
at  any  one  time  could  be  selected  from  a  library  of  such 
ROMs  by  the  Commanding  Officer  depending  on  the  ship's 
mission. 

The  determination  of  what  values  of  variance  of  accelera¬ 
tion  to  be  used  with  various  targets  are  best  determined  by 
modeling  Z~4_7.  A  computer  model  of  the  tracking  problem 
was  written  and  used  to  investigate  the  system  performance 
with  various  values  of  variance  of  acceleration.  The 
results  of  this  investigation  are  presented  in  the  figures 
discussed  in  Section  IV  E. 

C.  KALMAN  FILTER  INITIAL  CONDITIONS 

As  illustrated  by  Figure  5»  there  are  two  basic  loops  in 
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the  Kalman  filter  algorithm.  They  represent  the  recursive 
calculations  of  the  optimal  estimates  of  the  covariance  of 
estimation  error  matrix  (P)  and  the  target  state  matrix  (Z). 
An  initial  condition  must  he  established  for  each  of  these 
based  on  a  priori  information  of  the  measurement  process  and 
the  targets  of  interest. 

The  initial  condition  of  the  target  state  matrix  is 
equal  to  the  expected  value  of  the  initial  state.  The 
expected  values  of  the  initial  position  elements,  Xq  and 
Yq,  are  the  first  position  measurements; 

Xq  =  R(0)  sin  9(0) 

Yq  =  R(0)  cos  9(0) 

The  expected  values  of  the  initial  velocity  elements,  X^ 
and  Yq,  are  zero.  Thus,  the  filtering  of  the  target  state 
does  not  begin  until  the  second  measurement. 

The  initial  condition  of  the  covariance  of  estimation 
error  matrix  is  equal  to  the  covariance  of  error  of  the 
initial  state  estimate. 

Po  =  P  ^^0  -  ^0-7  -  7o-7  " 

E^(Xq-Xq)2^  0  0  0 

0  eZ"(Xq)2^  0  0 

0  0  E^(Yq-Yq)2J7  0 

0  0  0 

The  non-zero  elements  in  this  matrix  are  functions  of 
the  initial  relative  velocities  which  might  be  encoimtered, 
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and  the  measurement  errors  of  the  initial  position.  Assuming 

a  uniform  probability  distribution  of  relative  velocities 

from  -50  yd/sec  to  +50  yd/sec  (+  46  m/sec),  and  an  initial 

range  of  50f000  yds  (46,000  m),  the  variance  of  initial 

2  2 

position  error  is  12,000  yds  (10,000  m  )  and  the  variance 
of  initial  velocity  error  is  830  yd^/sec^  (700  m^/sec^) 
for  both  the  X  and  Y  coordinates . 

The  primary  influence  of  the  initial  conditions  is 
the  time  required  for  the  filter  to  achieve  a  steady-state 
solution.  Large  values  lead  to  longer  settling  times  for 
nonmaneuvering  targets  where  small  values  lead  to  longer 
settling  times  for  maneuvering  targets. 

D.  SYSTEM  COMPUTER  MODEL 

In  order  to  evaluate  the  performance  of  the  Kalman  filter 
with  various  targets  and  parameters,  a  computer  model  of 
the  system  was  derived.  The  computer  model  included  the 
simulation  of  target  and  own  ship  motion,  the  radar  measure¬ 
ment  and  interface  process,  the  Kalman  filter,  and  the 
computation  of  the  filtered  target  course  and  speed. 
Additionally,  the  program  included  computation  of  measure¬ 
ment  errors,  target  track  errors,  and  associated  statistics. 

The  computer  model  simulated  the  target  and  own  ship 
motion  by  digitally  integrating  the  north/south  and  east/ 
west  velocity  components  using  an  integration  step  of  0.01 
sec  and  double  precision  arithmetic.  Target  and  own  ship 
course  changes  were  accomplished  by  a  second-order  nonlinear 
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control  system  model  with  the  ordered  course  as  the  input, 
a  limit  on  the  maximum  angular  rate,  and  the  actual  course 
steered  as  the  output.  The  control  system  model  parameters 
and  angular  rate  limit  were  adjusted  to  simulate  typical 
course  changes  for  the  targets  of  interest.  Speed  changes 
were  accomplished  by  a  first-order  linear  control  system 
model  with  the  ordered  speed  as  the  input  and  the  actual 
speed  as  the  output.  No  attempt  was  made  to  duplicate 
actual  ship  maneuvers ,  but  only  to  approximate  the  maneuvers 
typical  of  the  targets  of  interest. 

The  radar  and  interface  systems  were  modeled  by  first 
comparing  the  rotation  of  the  radar  antenna  to  the  bear¬ 
ing  to  each  of  the  targets.  Upon  detection  of  a  target, 
the  bearing  and  range  to  the  target  were  noted  along  with 
own  ship's  course  and  speed,  and  the  time  of  detection. 
Random  numbers  representing  the  Gaussian  distribution  of 
measurement  errors  were  added  to  the  bearing  and  range 
measurement.  The  measurements  were  then  truncated  to 
simulate  the  digitizing  process  of  the  interface  system. 

The  Kalman  filter  was  implemented  through  FORTRAN 
statements  similar  to  the  algebraic  expressions  listed  in 
Appendix  A.  The  only  modifications  were  those  related  to 
improving  the  efficiency  of  executing  the  computations. 

The  target's  position,  course,  and  speed  were  computed 
from  the  relative  north/south  and  east/west  positions 
and  velocities  provided  by  the  Kalman  filter,  and  the 
simulated  measurement  of  ov/n  ship's  course  and  speed. 
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The  computer  model  was  programmed  in  FORTRAN  IV,  com¬ 
piled  with  a  FORTRAN  H  Compiler,  and  executed  on  the  Naval 
Postgraduate  School's  IBM  360/67  computer. 

E.  KALMAN  FILTER  PERFORMANCE 

Using  the  computer  model  previously  described,  the  per¬ 
formance  of  the  Kalman  filter  was  evaluated  for  two  classes 
of  targets,  maneuvering  and  nonmaneuvering.  In  all  cases, 
own  ship  was  proceeding  on  a  steady  course  (000°)  and  at 
a  fixed  speed  of  20  knots  (37  km/h)  from  an  initial  position 
at  the  origin  of  the  cartesian  coordinate  system.  Eight 
targets  were  initially  positioned  and  proceeding  on  the 
courses  and  at  the  speeds  shown  in  Table  1 .  A  unique  set 
of  random  numbers  was  used  with  each  target  for  the  various 
computer  runs.  This  eliminated  the  influence  of  variations 
in  the  random  number  distributions  on  the  statistics  of 
the  Kalman  filter  output. 

The  following  system  parameters  were  used: 


RADAR  PARAMETERS 


Pulse  Repetition  Rate 
Range  Precision 


100.0  pps 

7.057^  yds 
(6.453  m) 

0.1258  deg 

17.0  rpm 


Bearing  Precision 


Antenna  Rotation  Rate 


DIGITAL  INTERFACE  PRECISION  PARAMETERS 


Bearing  LSB 


Range  LSB 


4.0  yds 
(3.6576  m) 

0.02197  deg 
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TARGET  N/S  POSITION  E/W  POSITION  COURSE  SPEED 

NO.  ( vards/meters )  ( yards/me ters )  (degrees)  ( knots /kms ) 
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TABLE  I  INITIAL  POSITION,  COURSE  AND  SPEED  FOR  SIMULATED  TARGETS 


Time  LSB 


0.1221  msec 


Own  Ship  Course  LSB  0.02197  deg 

Own  Ship  Speed  ISB  O.03125  knot 

(0.05788  km/h) 

p 

The  Kalman  filter  measurement  error  parameters  (cri  and 

2 

0^  )  were  selected  to  match  the  computed  measurement  errors 
resulting  from  the  radar  and  interface  parameters. 

cf^  =  4.870  X  10"^  rad^  (0.01599  deg^) 

=  51.14  yd^  (42.76  m^) 

The  Kalman  filter  initial  conditions  were  the  same  as  those 
discussed  in  Section  C  of  this  chapter.  The  variance  of 
target  acceleration  (cr  )  parameter  was  chosen  as  the  inde- 
pendent  variable. 

Both  scenarios  began  with  the  acquisition  process  for 
each  of  the  eight  targets,  which  involved  setting  the 
initial  conditions  of  the  Kalman  filter,  and  ran  for  5 
minutes  of  simulated  time.  For  the  nonmaneuvering  case, 
no  changes  were  made  throughout  the  run.  For  the  maneuver¬ 
ing  case,  each  of  the  targets  executed  a  45  deg  course 
change  beginning  120  sec  after  the  start  of  the  run. 

Targets  1,  4,  6  and  8  turned  right,  while  targets  2,  3> 

5  and  7  turned  left.  These  turns  were  executed  with  an 
angular  rate  of  3*0  deg/sec.  Figure  6  shows  an  example 
of  the  plot  developed  during  such  a  left  turn. 

The  actual  measurement  error  statistics,  based  on  680 
measurement  samples  from  the  computer  model,  for  the 
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nonmaneuvering  and  the  maneuvering  target  scenarios  are 
listed  in  Table  II.  There  was  good  correlation  between  the 
measured  variances  of  error  and  the  computed  variances . 
Additionally,  the  mean  errors  were  approximately  equal  to 
minus  one-half  the  values  of  the  least  significant  bit 
used  in  the  interface  simulation.  This  is  in  agreement 
with  the  theoretical  mean  error  generated  by  truncation. 

Each  scenario  was  repeated  for  seventeen  logarithmically 
distributed  values  of  the  Kalman  filter  variance  of  target 

o  li  o 

acceleration  parameter  from  0.0005  yd  /sec  (0.00042  m  / 
i},  2  /  4*  /  2  / 

sec  )  to  1  yd  /sec  (0.85  ^  /sec  ).  In  all  cases,  the 
Kalman  filter  achieved  a  steady  state  condition  in  60  sec, 
corresponding  to  17  target  position  measurements.  Steady 
state  is  defined  as  the  condition  where  the  Kalman  filter 
gain  matrix  (K^)  would  be  constant  for  targets  with  constant 
relative  position  and  constant  sample  intervals.  The 
bearing,  range,  course  and  speed  of  each  target  was  com¬ 
puted  from  the  Kalman  filter's  optimal  estimate  of  target 
state  and  compared  with  the  actual  target  bearing,  range, 

t 

course  and  speed.  Error  averages  were  computed  for  the 
time  interval  from  60  sec  to  300  sec.  Figures  7,  8,  9 
and  10  show  plots  of  the  average  root  mean  square  of  error 
for  the  target  bearings,  ranges,  courses  and  speeds, 
respectively,  as  a  function  of  the  variance  of  target 
acceleration  parameter.  It  should  be  noted  that  the 
error  averages  for  the  maneuvering  target  are  a  function  of 
the  elapsed  time  between  maneuvers  and  the  time  spent  in 
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TABLE  II  MEASUREMENT  ERROR  STATISTICS 


CO  -H  O  ^ 
^  U  U  Q) 
O:;  Cd  ^  Td 

fq 


Cs] 
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FIGURE  7  AVERAGE  ROOT  MEAN  SQUARE  OF  BEARING  ERROR  VS. 

VARIANCE  OF  TARGET  ACCELERATION  PARAMETER  (o' 
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FIGURE  8  AVERAGE  ROOT  MEAN  SQUARE  OF  RANGE  ERROR  VS. 

VARIANCE  OF  TARGET  ACCELERATION  PARAMETER  (o' 
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FIGURE  9  AVERAGE  ROOT  MEAN  SQUARE  OF  COURSE  ERROR  VS. 

VARIANCE  OF  TARGET  ACCELERATION  PARAMETER  (a; 
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FIGURE  10  AVERAGE  ROOT  MEAN  SQUARE  OF  SPEED  ERROR  VS. 

VARIANCE  OF  TARGET  ACCELERATION  PARAMETER  (cr 


maneuver.  Thus,  the  plots  for  the  maneuvering  target  only 
apply  for  the  particular  scenario  used. 

Several  conclusions  can  be  dravm  from  these  figures. 

For  nonmaneuvering  targets,  the  variance  of  acceleration 

p  oh, 

parameter  (cr  )  should  be  small  (<  .001  yd  /sec  (.0008 

3. 

2  / 

m  /sec  ))  if  smooth  course  and  speed  values  are  desired. 

For  maneuvering  targets,  there  is  an  optiaml  value  of 
O' ^  (».10  yd^/sec^  (.08  m^/sec^))  below  which  the  filter 

cL 

does  not  follow  the  target  maneuvers  well  and  above  which 
the  target  data  will  be  noisy  regardless  of  the  target 
maneuvers .  There  is  a  lower  limit  to  how  well  the  Kalman 
filter  will  track  a  maneuvering  target  even  with  an  optimal 
value  of  O'  .  Finally,  the  bearing  and  range  to  a  non- 

3, 

maneuvering  target  will  be  improved  by  the  Kalman  filter, 
as  shown  by  the  standard  deviation  of  measurement  error 
line  on  Figures  7  and  8,  over  a  wide  range  of  o*  values. 

3. 
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V.  MICROPROCESSOR  PROGRAM 


A  MCS-8  microcomputer  system,  using  INTEL’S  8080B 
microprocessor,  was  used  to  implement  the  Kalman  filter 
algorithm.  The  objective  of  implementing  the  algorithm 
was  to  determine  limits  to  the  time  required  to  accomplish 
the  Kalman  filter  computations  and  the  accuracy  of  the 
results . 

A.  PROGRAMMING  LANGUAGE 

INTEL  Corporation  has  developed  the  PL/M  compiler,  a 
derivative  of  PL/l,  for  use  with  their  8080  microprocessor. 
PL/M  is  a  block  structured  language  which  speeds  micro¬ 
processor  program  development  by  relieving  the  programmer 
of  the  tedious  task  of  memory  and  register  management 
It  facilitates  the  use  of  variable  names  while  still  permit¬ 
ting  the  programmer  the  freedom  to  reference  specific 
memory  locations.  The  basic  instruction  set  includes 
arithmetic  and  Boolean  algebra  operation  on  8-bit  and  l6-bit 
variables  and/or  constants.  Program  control  is  accomplished 
with  Procedure,  Do,  Go  To,  If,  Call  and  Return  statements. 
Logic  functions  such  as  comparison,  shift  and  rotation  are 
also  included.  Floating  point  and  transcendental  functions 
must  be  implemented  by  subroutines .  PL/M  permits  the 
program  control  provided  by  assembly  languages  while 
taking  a  much  more  readable  form. 
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B.  FLOATING  POINT  ARITHMETIC 


Floating  point  arithmetic  is  required  whenever  the  range 
of  numbers  represented  by  the  program  variables  is  large  or 
unpredictable  Such  is  the  case  with  the  Kalman 

filter  algorithm.  This  is  particularly  true  if  the  operator 
is  given  the  freedom  to  alter  the  filter  parameters .  A 
disadvantage  of  using  floating  point  arithmetic  is  the 
additional  execution  time  required. 

A  suitable  set  of  floating  point  arithmetic  procedures 
were  not  available.  The  procedures  listed  in  Appendix  B 
were  derived  from  several  existing  procedures  and  the 
general  guidelines  detailed  in  Reference  12. 

A  l6-bit  mantissa,  7-bit  exponent  and  1-bit  sign  were 
chosen  to  represent  the  floating  point  numbers  using  three 
8-bit  data  words  in  the  following  format: 


7  0  7  0  7 _  0 


M 

M  +  1 

1  M  +  2 

Mantissa  (l6-bits)  '  Sign  (1-bit) 

Exponent  ( 7 -bits ) 

M  -  Memory  Address  Associated  with  the  Floating  Point 
Number 

The  l6-bit  mantissa  was  chosen  because  it  is  the  largest 
binary  nimiber  that  INTEL'S  8080  microprocessor  can  manipu¬ 
late  without  a  considerable  increase  in  the  comjTLexity  of 
the  progi^Linming .  Additionally,  it  is  compatible  with  the 
input  and  output  data.  The  mantissa  is  always  left  justi¬ 
fied  with  the  binary  point  located  to  the  left  of  the  most 
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significant  bit.  Thus  the  mantissa  is  always  equal  to  or 
greater  than  0.5»  and  less  than  1.0. 

The  sign  was  located  in  the  third  word  in  order  that  the 
mantissa  be  a  full  l6-bits.  A  positive  sign  is  represented 
by  a  binary  0  and  a  negative  sign  is  represented  by  a 
binary  1 . 

The  remaining  7-bits  of  the  third  word  are  used  for  the 
exponent.  The  value  of  the  exponent  is  equal  to  64  plus  the 
power  of  two  to  be  associated  with  the  mantissa. 

The  data  and  constants  in  the  PL/M  programs  are  written 
in  this  format  with  hexadecimal  numbers.  The  following 
are  examples  of  this  form. 


DECIMAL 

0.5 

1.0 

-1.0 

10.0 

64.0 

0.1 


HEXADECIMAL  FLOATING  POINT 

80H,  OOH,  40H 
80H,  OOH,  41H  ^ 

80H,  OOH,  CIH 
AOH,  OOH,  44H 
80H,  OOH,  47H 
CCH,  CDH,  3DH 


The  largest  magnitude  number  that  can  be  represented 

is  FFH,  FFH,  7FH,  which  corresponds  to  9-223  x  10^®.  The 

smallest  magnitude  number  that  can  be  represented  is  80H, 

—  20 

OOH,  OOH  which  corresponds  to  2.710  x  10  .  Zero  is 

represented  by  OOH,  OOH,  OOH. 

The  typical  execution  time  of  the  arithmetic  and  trig¬ 
onometric  procedures  listed  in  Appendix  B,  using  INTEL'S 
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8080B  microprocessor,  are  listed  below. 


FLOATING  POINT  PROCEDURE  MAXIMUM  EXECUTION  TIME  (msec) 


Addition  (ADD)  1,3 
Subtraction  (SUB)  1.4 
Multiplication  (MULT)  3.4 
Division  (DIV)  4.1 
Comparison  (COB/CPARE)  0.8 
Square  Root  (SQRT)  17*0 
Cos  (TRIG  1)  18.0 
Sin  (TRIG  2)  18.0 
Cos  and  Sin  (TRIG  3)  30.0 
Arctangent  (TRIG  4)  19.0 


The  actual  execution  times  of  these  procedures  are  a  fionc- 
tion  of  the  numbers  involved. 

Overflow  and  underflow  provisions  are  included  in  the 
arithmetic  procedures.  Overflow  will  result  in  a  solution 
equal  to  tlie  largest  possible  number  with  the  appropriate 
sign.  Underflow  will  result  in  a  solution  equal  to  zero. 

l/jC'jiHi  tr 

The  precision  of  the  sine  and  cosine  functions,  using 
the  procedure  "TRIG”,  is  at  worst  +  0.000488  if  the  argu¬ 
ment  is  in  the  range  0  to  Str.  The  precision  of  the  arc¬ 
tangent  function,  using  the  procedure  "TRIG",  is  equal 
to  +  0.0610  mrad. 

The  four  arithmetic  procedures  are  used  by  passing  the 
variables'  addresses  to  the  procedure  with  a  call  statement. 
The  following  PL/M  statements  are  examples. 
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FUNCTION 


pl/m  statements 


X  +  Y  =  Z 
X  -  Y  =  Z 
A  X  B  =  C 
L  f  M  =  Q 


Call  ADD  ( .X,  .Y,  .Z) ; 
Call  SUB  ( .X,  .Y,  .Z) ; 
Call  MULT  ( .A,  .B,  .C) ; 
Call  DIV  ( .L,  .M,  .Q) ; 


Note:  The  dot  symbol  proceeding  a  variable  name  in  PL/M 


indicates  the  address  of  the  variable.  Some 
variables  may  correspond  to  addresses  of  other 
variables,  in  which  case  the  dot  would  be  omitted. 


The  "COMPARE"  procedure  returns  a  binary  integer  quantity 
which  indicates  the  relationship  of  the  arguments . 

pl/m  statement 

N  =  COMPARE  ( .A,  .B) ; 

RESULT 

If  A  <  B  Then  N  =  0 

If  A  =  B  Then  N  =  1 

If  A  >  B  Then  N  =  2 


The  argument  for  "SQRT"  must  be  positive. 

FUNCTION  PL/M  STATEMENT 

Y  Call  SQRT  ( .X,  .Y) ; 

The  procedure  "TRIG"  can  be  used  to  compute  the  sine 
and/or  cosine,  and  arctangent  by  passing  the  addresses  of 
the  arguments  plus  a  function  indicator  flag  to  the  pro¬ 
cedure  with  a  call  statement.  It  requires  less  time  to 
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compute  sine  and  cosine  of  the  same  angle  with  one  call 
rather  than  two.  The  following  are  examples  of  the  multiple 
uses  of  the  "TRIG"  procedure. 


FUNCTION 
X  =  Cos  (A) 

Y  =  Sin  (A) 

X  =  Cos  (A)  and 

Y  =  Sin  (A) 

A  =  Arc tan  (Y/X) 


PL/M  STATEMENTS 

Call  TRIG  ( .X,  .W,  .A,  1) ; 

Call  TRIG  (.W,  .Y,  .A,  2); 

Call  TRIG  (.X,  .Y,  .A,  3)j 
Call  TRIG  (.X,  .Y,  .A,  4); 


Note:  W  is  a  dummy  variable  used  to  fill  the  unused  argu¬ 

ment  location  in  the  call  statement  when  only  cosine 
or  sine  is  desired. 


C.  TRACKING  SYSTEM  PROGRAM 

An  experimental  version  of  the  tracking  system  software 
was  written  to  determine  the  execution  time  of  the  Kalman 
filter  algorithm  and  the  error  statistics  of  the  resulting 
target  track  data.  This  program,  written  in  PL/M,  is 
listed  in  Appendix  C.  There  are  three  functional  parts 
to  this  program;  (1)  data  input,  (2)  the  Kalman  filter,  and 
(3)  computation  of  output  data. 

The  data  input  portion  of  the  program  includes  the 
simulated  interface  system  output  data,  and  the  simulated 
control  panel  functions.  The  simulated  interface  system 
data  was  generated  in  the  correct  floating  point  format 
by  the  FORTRAN  program  used  to  model  the  system  performance. 
The  variables  "BRG,"  "RNG, "  and  "TIM"  correspond  to  the 
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measured  target  bearing,  range  and  time  of  detection, 
respectively.  The  ovm  ship  north/south  and  east/west 
velocity  are  represented  by  the  constants  "NSVEL"  and 
"EWVEL.  "  The  Kalman  filter  parameters  cr  ^ and  cr  ^ 

SL  w 

are  represented  by  "SIGA,"  "SIGTH,"  and  "SIGR."  The 
variables  "MARKOLD"  and  "CONTROL"  serve  the  functions  of 
distinguishing  which  target  data  corresponds  to  existing 
tracks,  and  which  control  panel  switches  are  turned  on. 
These  variables  are  associated  with  the  variables  "SIMFLAG 
(  ) "  and  "DFLAG"  which  correspond  to  the  interface  unit 
from  which  the  target  data  came  and  thus  the  target's  track 
number.  In  an  actual  operating  system  the  program  would 
use  INPUT  statements  or  direct  memory  access  to  read  in  the 
target  measurement  and  control  panel  switch  data.  The 
measurement  data  would  have  to  be  converted  from  fixed 
point  to  floating  point  format  before  processing. 

The  Kalman  filter  statements  are  grouped  into  two  long 
DO  statements.  The  first  is  used  to  initialize  the  filter 
when  a  new  target  is  first  processed.  The  second  group  is 
used  to  process  new  data  for  existing  target  tracks.  The 
sequence  of  statements  and  variable  names  parallel  the 
algebraic  expressions  listed  in  Appendix  A.  Temporary 
variables  were  introduced  where  necessary  to  avoid  repeti¬ 
tive  operations  and  to  store  intermediate  results. 

Using  the  Kalman  filter  optimal  estimates  of  target 
position  and  relative  velocity,  the  third  portion  of  the 
program  computes  and  displays  the  target  bearing,  range, 
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course  and  speed. 

Following  the  program  debugging,  a  series  of  tests, 
using  a  full  range  of  typical  input  data,  were  run  to 
verify  the  correct  computation  of  each  variable.  The 
FORTRAN  model  was  used  to  compute  the  "correct"  values  for 
each  variable  using  the  Naval  Postgraduate  School  IBM  360 
computer.  It  was  found  that  the  computational  errors  were 
confined  to  a  range  of  significance  sufficiently  small  to 
avoid  instabilities  in  the  filter.  After  processing  several 
measurements,  the  Kalman  filter  estimates  of  target  state 
computed  with  the  microprocessor  were  within  the  same 
range  of  error  as  the  IBM  360  solutions.  The  covariance 
of  estimation  error  and  Kalman  gain  matrix  elements  computed 
by  the  microprocessor  converged  rapidly  on  the  values 
computed  by  the  IBM  3^0  computer. 

As  a  check  on  the  microprocessor's  performance  a  series 
of  eight  test  runs  were  made,  each  with  a  different  target. 
The  eight  simulated  targets  were  identical  to  the  non¬ 
maneuvering  targets  used  previously  except  the  signs  of 
the  east/west  initial  positions  for  targets  3-8  were 
negative.  This  change  was  made  because  the  execution  time 
and  precision  of  the  "TRIG"  procedure  is  worst  for  targets 
in  the  north-west  quadrant.  Nonmaneuvering  targets  were 
used  because  the  available  memory  space  limited  the  amount 
of  data  that  could  be  loaded  with  the  program.  The  test 
runs  were  restricted  to  I30  sec  tracks.  A  variance  of 
acceleration  of  .001  yd  /see  (.00084  m  /sec  )  was  selected 
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to  insure  that  computational  errors  would  not  be  over¬ 
shadowed  by  normal  estimation  error  in  the  output.  It 
was  assxamed  that  these  tests  would  represent  worst  case 
conditions . 

The  test  program  occupied  a  total  of  2E00^^  bytes  of 
memory.  This  was  divided  as  follows; 


Floating  Point  Procedures 

98^16 

Bytes 

Kalman  Filter  Algorithm 

12C7ig 

Bytes 

Bearing,  Range,  Course  and 
Speed  Computations 

22O16 

Bytes 

Variable  Storage 

383i6 

Bytes 

The  error  averages  for  the  test  data  are  listed  in 
Table  III.  These  error  averages  represent  the  performance 
of  the  microprocessor  tracking  program  after  the  Kalman 
filter  achieved  a  steady  state  condition,  approximately 
50  sec  after  initialization.  This  resulted  in  a  relatively 
small  sample  set  of  20-21  samples  for  each  target.  The 
error  statistics  of  the  input  data  during  this  time  are 
listed  in  Table  IV.  As  compared  to  the  measurement  error 
statistics  presented  in  Table  II,  there  is  less  correlation 
with  the  theoretical  measurement  errors  for  the  microprocessor 
test  data.  Therefore,  the  error  averages  listed  in  Table 
III  should  be  considered  typical  only  for  a  short  period 
of  operation  rather  than  for  long  periods  with  all  possible 
targets . 

An  analysis  of  these  results  reveals  the  microprocessor's 
influence  on  the  track  data.  By  comparing  the  course  and 
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TARGET  BEARING  RANGE  COURSE  SPEED 

NO.  (degrees)  (yards)  (degrees)  (knots) 

_  Mean _ RMS  Mean _ RMS  Mean _  RMS  Mean  RMS 
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TABLE  III  MICROPROCESSOR  TRACKING  PROGRAM  ERROR  AVERAGES 


TARGET  BEARING  RANGE 

NO.  (degrees)  (yards) 

Mean  Std.  Dev.  Mean _ Std.  Dev 
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TRACKING  PROGRAM  INPUT  DATA 


speed  data  for  targets  3i  5  and  7  with  the  data  for  targets 
4,  6  and  8  it  can  be  seen  that  for  targets  with  similar 
course  and  speed  the  errors  increase  with  range.  Addition¬ 
ally,  errors  are  larger  for  targets  with  small  relative 
velocities.  The  combination  of  long  range  and  small  rela¬ 
tive  velocities  can  be  expected  to  result  in  larger  errors 
due  to  the  limitations  of  the  floating  point  arithmetic  to 
handle  large  numbers  with  small  differences. 

The  final  parameter  of  interest  is  the  execution  time 
of  the  microprocessor  program.  The  INTEL  8080B  micro¬ 
processor  was  able  to  execute  85  consecutive  iterations  of 
the  Kalman  filter  algorithm  in  37*5  secs,  resulting  in  an 
average  execution  time  of  0.44  secs.  The  computation  of 
target  course  and  speed  required  a  worst  case  execution 
time  of  54  msec.  It  would  not  be  necessary  to  compute 
the  target  course  and  speed  for  every  target  on  each  antenna 
scan.  The  additional  processing  required  to  output  the 
course  and  speed  data  would  best  be  accomplished  by  a 
separate  microprocessor  dedicated  to  each  display  device. 

In  general,  a  network  of  microprocessors  could  be  used  to 
share  the  task  and  thus  increase  the  number  of  targets 
which  could  be  tracked. 
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VI.  CONCLUSION 


The  tedious,  error-prone  task  of  tracking  surface  search 
radar  targets  manually,  as  is  common  onboard  naval  ships 
not  equipped  with  NTDS,  can  be  eliminated  with  the  use  of 
large  scale  integrated  microprocessors. 

The  AN/SPS-10  Surface  Search  Radar,  installed  on  most 
U.S.  Navy  ships,  could  provide  target  position  measurements 
with  sufficient  accuracy  to  perform  automatic  tracking  of 
surface  radar  targets.  Using  leading  edge  detection  and 
currently  available  hardware  it  would  be  possible  to 
digitally  measure  the  bearing  and  range  to  a  target  with 
a  standard  deviation  of  error  of  0.126  deg  and  7*15  yards 
(6.5^  m)  respectively. 

Using  a  4-state  Kalman  filter  and  appropriate  conver¬ 
sions  for  polar  coordinate  measurements,  it  is  possible  to 
obtain  optimal  estimates  of  target  north/south  and  east/ 
west  position  and  velocity.  From  this  data  and  knowledge 
of  own  ship's  motion,  target  course  and  speed  can  be  com¬ 
puted  far  more  accurately  than  by  manual  tracking  tech¬ 
niques  . 

A  set  of  floating  point  arithmetic  procedures  was 
developed  for  INTEL'S  8080  microprocessor  based  on  a 
floating  point  number  format  with  a  l6-bit  mantissa  and 
a  7-bit  exponent.  These  procedures  maintain  sufficient 
accuracy  to  reliably  compute  the  Kalman  filter  variables 


79 


while  limiting  the  execution  time  sufficiently  to  pennit 
complex  functions,  like  the  Kalman  filter,  to  he  implemented. 

The  Kalman  filter,  and  course  and  speed  calculations  were 
successfully  implemented  using  the  microprocessor  develop¬ 
ment  system  at  the  Naval  Postgraduate  School.  ^  The  average 
execution  time  of  the  Kalman  filter  was  0.44  secs  while 
the  course  and  speed  calculations  required  54  msec.  This 
was  accomplished  with  computed  course  and  speed  errors  for 
nonmaneuvering  targets  typically  less  than  2.0  deg  and  0.5 
knots.  With  the  3*5  sec  sweep  of  the  AN/SPS-10,  a  single 
system  could  continuously  provide  course  and  speed  data 
for  seven  targets.  This  could  he  increased  hy  using  faster 
microprocessors  currently  entering  the  market. 

It  is  recommended  that  future  work  in  this  area  he 
devoted  to  (1)  improvements  in  the  floating-point  software, 
(2)  development  of  the  interface  hardware,  and  (3)  improve¬ 
ments  to  the  Kalman  filter  algorithm.  The  floating  point 
software  could  he  improved  hy  reducing  the  execution  time 
of  the  arithmetic  procedures  and  increasing  the  precision 
of  the  transcendental  functions.  The  Kalman  filter  algorithm 
could  he  improved  hy  minimizing  the  effects  of  own  ship's 
motion,  and  reducing  the  execution  time. 

A  microprocessor-hased  tracking  system  has  the  potential 
to  provide  any  ship  with  a  digital  automatic  surface  track¬ 
ing  capability,  without  the  expense  of  NTDS  size  equipment. 

If  implemented,  it  could  reduce  the  manning  of  underway 
watches  while  improving  the  quality  of  the  tactical 
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information  available  to  the  Officer  of  the  Deck.  There 
are  additional  applications  of  similar  systems  for  such 
functions  as  sonar  tracking  and  navigation. 
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P(n)  =  i  P(n-1)«^  +nQr 
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APPENDIX  A  ALGEBRAIC  FORM  OF  KALMAN  FILTER 
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K(n)  =  P(n)  H^(H  P(n)  +  R(n)  )“ 
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P(n)  =  (I  -  K(n)  H)  P(n) 
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APPENDIX  B  FLOATING  POINT  SOFTWARE 


M 

<l 


^ 

lU 

iU  U1 

o  z 


2 

O 

•* 

•  * 

<l 

x-4 

U 

i/) 

H- 

h- 

f-H 

•«k 

u. 

^L 

X 

X 

»»4 

U4 

UJ 

UJ 

■r 

•  «k 

X 

— 

3 

h- 

i-J 

a 

fVJ 

1— 

1— 

H- 

w 

>v 

1^ 

CO 

— • 

X 

-M* 

X 

X 

O 

O 

It 

X 

LL. 

o 

o 

o 

— 

> 

O 

GO 

og 

o 

:z 

X 

X 

1- 

UJ 

- o 

lU 

UJ 

o 

X 

X 

o 

o 

h- 

K 

Q  ••» 

►- 

o 

<l 

x>- 

o 

X 

CL 

3 

3 

Xr-4 

II 

UJ 

-•X 

CtCCCvJO 

•  •  <1 

•  •  <1  •» 

UJ^ 

— * 

"Z. 

X 

2 

Od'—  H 

O  3f\J 

a  om 

XX 

fVJ 

o 

XX 

•-< 

Q  CJH- 

Q  C/H- 

1—  It 

X 

X 

— • 

'Nl 

X 

X 

1- 

V.COCO  II  3 

* 

3— tox 

X-*coX 

3*-l 

X 

<13 

1— < 

.. 

•  •* 

«xcvj2:uj 

UJCVJXUJ 

•  » 

•  • 

x<  •• 

3 

XUJ 

-JO 

XwO 

XwO 

OisJX 

o 

AO 

••3fO 

i/lX>OcO 

<IQ 

h-x»-«o 

O  QC 

It 

Q 

>-Oh- 

UJ 

<XwwC0Z 

3 

u?h- 

II  COH- 

X 

>-X»-» 

cc 

COUJUJ  o 

X  tl 

X 

coxi— 

X  - 

*  LO  A 

<l 

CO  J  JO»-» 

OIUJ 

XaiLLO 

XUJLUO 

X  II  UJX'-*«*‘ 

XUJ^ 

XX 

♦-•cDcaz^ 

X 

VMhhO 

Of^QCUjrsI^ 

0X04X0 

i/I 

i<i 

i/ll— 

>- 

cx. 

f~W 

’X«3 

Ul 

zao  UJ 

1- 

> 

>•  •• 

cO 

•^3  •• 

h-  Xi/I'/JK-** 

-1 

-•aacj^co 

ZLU 

a 

Q 

0»-a 

X*^  tl 

It  O 

CD 

LUX 

>s.LL  X 

V.X  X 

It  UJX 

ai>-iu>-xo2: 

<3. 

^  tl 

•rt*  *-4  UJ 

•1^  H-t  LU 

rvJUC.3 

XA»w>-»-houj 

^xxx 

ljlu 

3UJ 

UC 

uuc/ivo*^ 

U.  A 

<l 

O  tl  II 

X 

A 

A 

II 

X 

a.  >0-0. 

uux 

X 

X 

UJX 

r»» 

•ft* 

•R- 

•R* 

's. 

89 


•ft-  -«• 


LU 

X 

h- 

•  «k 

»u 

X 

X 

<1 

fNl 

M 

rsj 

•  «k 

11 

11 

II 

•ft- 

•ft- 

H- 

UJ 

IM 

IM 

«M 

w 

X 

> 

* 

X 

UJ 

M 

A 

M 

A 

PsJ 

lU 

oc. 

>- 

X 

X 

<1 

•  «k 

«•» 

•  «k 

o 

O 

.  Q 

LL 

h- 

CO 

fN| 

z 

•Nl 

z 

1X1 

CD 

< 

IV| 

•«x 

PSi 

< 

INI 

o 

w 

w 

H* 

o 

3: 

1— 

3 

h- 

3 

Z 

►-* 

o 

z 

a 

z 

O 

UJ 

z 

H- 

-j 

UJ 

X 

X 

z 

••*UJ 

z 

II 

OL 

II 

CL 

II 

a 

-'X 

LU 

X 

a. 

OJH- 

as.  ^ 

t-i 

X 

f-4 

X 

X 

U) 

w 

X 

X 

UJ 

Z-J 

lUQ 

PsI 

1— « 

rsl 

1—4 

PnI 

3CUJ 

Q 

Q 

•  •k 

A 

a:3cs4 

j-Q  O— 

•  •» 

1— 

CD  CJ«-I 

— • 

X 

UI-— 

X 

CO^ 

•  XUJftH 

LL 

rsj 

QC. 

••-OPvl 

a: 

•-Xrxl 

(MXc/>X 

O  Z»^*-< 

•Nj 

<l 

X-DPsI 

>-“^rxl 

U.W 

’^ZUJ 

LU  LUrsI  -f 

XCD^  — 

o>cx:o 

lOOC  X-* 

xz 

CO 

1  cxz 

</) 

1  <ixz 

X»-»UJ 

1— OC. 

CDOL 

z 

^  ox 

z 

X  ox 

H- 

CD  OUJ 

o 

>-X^-D 

e> 

XX4-»ID 

z  n 

II 

I-4CO  >-.^COMaX(— 

1—4 

II  XXI— 

1-^ 

II  XXH- 

UJLUXLLO 

vjQc.  >cr  II  II  z  II  ui 

PvJ<t  II  X 

CO 

M<1  11  X 

UJ  cx^^UJUJiNirx 

rsiorslOL 

iNH_;rNIX 

ZOi  ^^FNKXI 

z 

z 

UL. 

UJ^ 

X 

X 

•  •k 

X 

JL^  II 

X 

X 

(-1 

UJ 

:«Z**'NJLL 

JC  •• 

x«* 

X 

IM 

■ft 

■ft  -ft  r— 

•ft-  f— 

•ft- 1- 

V^I-4 

CD 

X 

X 

X 

z 

lU 

X 

X 

X 

90 


FLC^IING  POINT  SUBTRACTION  ROUTINE 


I/) 

CO 

tu 

a: 

•♦Q 


M  >>- 
•m>^(DCD 

<<t 

>-f>4<  — 
•»  •‘>•00 
<<1  w 
X>-CJ 
W  •‘OJCO 

<IC03 

Z>  X 
Qai>-> 
UJQL^^ 

o<*. 

LJ— I 
CXO 
U.UJ 
CJ 


X— 

o< 

COfNl 

xco 

03 

X2 


*-IC\J>- 


>->•  • 
>->-  •• 

M  II  < 

II  w«i«rCJ 

^yicocoo 

i.5S-X-JUJ 

>->->-Ot/> 

CJ 

2: 

OJ 


FLC^TING  POINT  MULTIPLY  ROUTINE 


X 

o 

<r 

I 


X 

o 

GU 

LJ 

<l 


Q 

X 

2? 

X 

L2J 

f^ 

CM 

•  •k 

Q 

w 

X 

>- 

ac 

<1. 

-J 

CC 

o 

•  * 

UJ 

CM 

X 

CO 

QC 

w 

CO 

>- 

UJ 

9m 

•  •k 

w 

CM 

oc 

o 

X 

•  •k 

w 

Q 

\  II 

X 

0 

+ 

X 

O 

* 

4* 

••  II 

w 

.•CV4 

•-4f— 1 

PNJ 

— 0-- 

w 

ao'-^ 

ww 

rNi 

XOCM 

QCOfNi 

x>- 

II 

CL 

or 

>- 

LU 

rsj 

r>-XM 

o 

fNl^  •» 

cx:cc 

rsi 

UJ 

— 1»  •‘IJJ 

cuo 

CIO 

ax 

UJ 

<^Xh-aiujLij 

oox  II 

Z 

21- 

^ - 

••UJ 

•  » 

^  ••  •• 

•» 

f-H 

GOOD 

O** 

^  xooto  — 

<<  CCCDCQ 

a:--.-* 

o  •* 

,  1.41^^ 

Q-- 

— U.  II  2 

•  •f*s| 

OIOM 

•k  f>-4 

f— 4  •• 

cMr— •111  II 

•*  pi>  <I  ^ 

oi  II 

UJX>- 

2r  •* 

»k^ 

rM'—' 

PNl'.^w 

UJINJ  1 

COXA'^ 

— rsj 

X>-UJ 

Z^^KD 

1  UJLU 

'-•iXrNi 

il 

1— w  M 

CM— ^ 

w 

2:  11 

2:-j-j 

'J-qcq:^ 

l-w  iU 

2  UJ  •• 

w  t3  •• 

UJ<C<^UJUJLU 

cr:r^ 

OGQCQ 

.-40<q: 

CC  M 

UJ  II  >  2rMX-J2 

QLXc/iVjviv; 

QCLJ 

•• 

2  <5.  CL 

c;  II  cx 

UJ 

OLJ** 

U  wCO 

•  » 

>-COO**UJ 

CJtUin  ••O 

CD  CDCOCO 

X-^ 

0?QQX 

1—  II  II 

X 

QC  II  II  0(0 

Q->  OK- 

X-IK- 

UJUJOt 

i-O 

QLw^U- 

>-a-rNl 

X 

IX^^CUZT-J 

X<X  II  UJ 

II  — UJ 

»-•  II 

UJ-J-JU. 

-h 

<].r>jrNiujiij 

ujco*-tf^ac 

pNifNicr. 

LUX 

OJXXU. 

II 

fNI 

UJ 

ix.^^ 

•-M 

nj 

UJ 

1 

U.CJO 

U. 

^  li  II  II 

UJ 

it 

UJ 

UJUJ 

ZX>-fNl 

CJ  2 

rsi 

u. 

<l 

•  •  t_jUJ 

LJ  VU 

rk^ 

92 


FLCATING  POINT  DIVICE  ROUINE 


CO 

90> 

CO 

X 

h- 

UJ 

0 

cr 

00 

X 

CJ 

X 

0 

0 

<i 

0 

<l 

h— 

A 

X 

a. 

0 

0 

X 

II 

CM  — 

0 

UJ 

— 

— * 

H- 

CM 

>- 

•• 

w 

^  x> 

•  «> 

m 

>- 

*-4 

> 

>- 

CL  CCCL 

•h 

>• 

0  CJO 

X 

— xaiujuj 

X 

1^ 

QC 

CXh-f-H- 

X 

—  — 

11 

X 

fNi  •» 

0 

—  OOCD 

X 

CO  •• 

good  CD  UJ 

II  ••• 

X 

•  •ktsj 

11  ^ 

<M  h- 

c 

>  + 

>  *•<<!<>-••* 

ODXX  X>- 

Q 

>x 

•^X>-MCDUJ 

0— ^UJ 

XX  ww 

•«k 

INI 

^>-  K- 

IVI 

^X 

2:  11 

X  ••CJCJCJ— > 

2:  •• 

Z  X_l-J 

xo 

XX 

-^<UJUJUJ2^CD 

UJ  21 

-QQQOCDCO 

xo 

X^sl 

X  ••x 

1— w 

f— 

QC^<<*.<CiO— 

l-XZ) 

— ^<  -^oo 

X 

:d  ojojoj  •<»> 

Oh- 

CM  XOO 

>•  11  Q 

xco 

QUJ  »-l^ 

X  II  UJ 

w**  <— » 

>-fNl2^ 

XX 

UJQC 

OrMCC 

-?^IMIM  1  -J— 1 

11  P»JX 

/NX 

II 

w^-^UJXX 

LJ^ 

11 

UCO 

^  11  II  11  tl  If 

a.uj 

X 

oxxxx>- 

X 

X 

LJ 

1-^ 

i-4 

> 


X 

cox 

II  X 


X 

2:  CO 
UJ  11  •• 
•••XX  — 
CJI— • 
a  •• 
XLUfsi 
ZTOCOINI 


UJOO_l^ 


XAXX  XX 

> 

1— X 

X  XX 

>- 

X 

co*^  coco  •• 

•  «.  1 

II  X  ••^  II  II  X 

>ox 

OCX 

iNi2:ofNix2: 

—IX 

II 

QC—I 

<l 

^slXQ^slXX 

OQ- 

1— X- 

X 

X  »«k 

CO  Q 

X 

X 

X  2 

X  X 

X  11 
II 


^^JU 

r>iLJ 


C-i 


93 


=  SHL(Z2,  U  +  l 
=  ShL(TEf<F,l  ) 


OVEFFLQW/UNCERFLOW? 


X 


a: 

O 

2: 

o 

to 

II 


rvi 


rvj 


O 

il 


X  CM 


U. 

li. 

o 

II 


.•.QM 

a 

QZ 
UJ  •«» 

zxx 

UJH-U. 


o 

II 


IXI 


1*  IX.  •• 

^-ujO(-J**oo 
>-  II  2:0  11  z  •• 
XAiNlUJOi^tUZ 
U-UJ  (X 


r-x 

LU 

A 

CO 

1-0 

LUO. 

ujz 

UJ 

ULUJ 

PM 

M 


u 

—I 

II 


rsiLU 

f>4 

••Ct 
— UJ 

MZ 


O  II 


X 

LU 


XlM 
II  '-'liJH-l 
r>4*>4a:.UJ 


t-J 

Z 

UJ 


TIFN 


FLC/STING  POINT  CQNPAPISGN  ROUTINE 


(NJOCMO 


•» 

zzzz 

</> 

x>- 

QCQCcca: 

••*(/) 

3X33 

UJUJ 

Qco: 

J— J— h-f— 

Kcr 

00 

LUlUtULU 

>0 

•  «k 

QCOCQCQC 

zzzz 

CM 

OIQ 

0 - 

••k 

acococac 

Oi^CVJ 

<  'JJ 

0  GOOD 

ozzzz 

3333 

z 

II  II  II 

—  h- 

•>  •» 

CJ  LULU  LULU 

QC 

LULULU 

<1—*  •• 

xxxx 

lUUJUJUJ 

3 

(racer 

>>IJUtUCQ 

UJXXX>- 

Zl~‘l"*h"l— 

QCLrOCUC 

1— 

«< 

•« 

XXX^'— 

m 

X 

0,0.0. 

<  •.>•>—* 

XX 

»— f^r^oiuj 

'^ZZZZ 

QC 

xxs 

XXCOcOZ 

00 

00 

rr  tULUXiUi 

X 

ooo 

0 

coco 

— OQcDcO 

AVAVZ 

xxxx 

ooo 

•»<!<(  I-H 

>vZZZx:d 

lUOUJiUXXoC 

XI— »— 1— 1— 

X 

UJ<X^CO 

OQ 

•»o<<oo 

a>  II  XXXX3 

OX 

cr>  ^ 

zz 

•-M  OQ 

l-HZ  H. 

<ll- 

V  II  /N 


a<cujuu2: 

o  aacut/) 
rrai  X 


OOsjCNU^J 
II  J- 


t— OU-XXXLU 

»i  M  QC 

VIVJ 


<tA  V  A  VZ 

oujujxxac. 


LUO 

CX.LU 

X 

X 

o 

o 


la 

oo 

tu  II  • 


IMCM 

LJZXX^IL/^ 

ox 

XXXXX3  ••* 

zzo 

*—»*■** 

XO  It  II  It  It 

X 

Z  1-0 

00 

x>- 

hhXXX>- 

X 

z 

xxxxxz 

i-iZx 

It  M 

vj(/;x^x>- 

L/;»— t 

X 

M  ■<i-Hcr  X 

xi/)acx 

zz 

z>^ 

z 

z 

zxo^ 

00 

0 

0 

0 

I-4XXX 

LOLO 

vix 

(/> 

l/J 

i^ii-HLro 

0 

'll* 

*«■ 

\  CJ 

95 


FLC^TING  POINT  SQU/»REROGT  ROUTINE 


cc. 

OJ 

CD 


If  -fr  ^ 


<I 

UJ 

CL 

UJ 

> 


00 

o 

CL 

OO 


UJ 

-J 

cu 

<r 


QC 

> 

UJ 

2: 


U1 

00 

<1. 


•R*  if  -Jf 


If  -ft* 


i/) 

•  mt 

►-4 

X 

0 

f— 

'T 

a 

0 

— 

cc 

•«» 

•» 

LL 

UJ 

0 

fsl  •• 

w  — «• 

•  •k 

•k 

CsJ 

X  04 

CO 

cc^ 

OJ 

CQI— 

CO 

CJUL 

X  •* 

CC 

•  • 

Ul 

^-4X 

c/;ai 

UX-« 

•k  •k 

•  onr  iJ  iLLi 

H- UJ 

If  OJ 

OCJ 

1— f- 

— OKI- 

<r 

—  1 

— O-J 

•  • 

<Q>-> 

OJw 

000  n 

•k  «k 

ivKcOCO 

< 

wq: 

•kOO-- 

<((]Q  ••k 

•k 

hhCsJ 

h-X 

X 

•* 

<r-— 

00 

xcc— 

1 

k 

X<C?<ixj  UJUJ 

0-- 

2  1 

wUI- 

— • 

wfNI 

QCh- 

UJX 

lU 

j>arj 

'k— 

Q-^ 

xo 

j— 

h* 

UJ<IUJUJUJGUCU 

J. 

Q.< 

f— *4“ 

aj»-4*« 

II 

catxv/v;f- 

9 

<3.^ 

M 

II 

<r 

X  — 

004  og 

LJ-J-J— •• 

cj  GucDoarnm 

1 

Orvj 

QoslOkJ 

H-— f-jrgu 

^<M 

LMliJ  wv 

co^ 

W'kn^'kn^ 

IN|w 

0->  Q1 fNl  €X.  UJ  f— 

H-# 

Vf- 

-JctX 

UJ 

f— 

h- 

UJ 

XXCII 

It 

II 

UC— « 

r^uj 

tJU 

•  •k^ 

CLW 

li 

U} 

M  II  X 

tu 

UJ 

U - 1 

04M  II 

CJ 

II  ^ 

••  LJ 

:«•  -ff 

r— <uj 

UJ 

O. 

to 


t  f 
— 

h-OL 

UJO^ 

OLwri 


96 


END 


FLC^TIKG  PGINT  COSINE  FUNCTION 


*»•**-»«•* 


CM 


Q. 


t| 

V 


< 

UJ 


II 

V 

o 

CJ 


X 

►  ► 

0 

XX 

► 

^4" 

00 

X  •» 

•‘X 

••'•T'T 

ox 

^xx 

X 

XOr«i 

XXX 

X  •‘'T 

00 

OILJCM 

u-x  •* 

^r  t-i 

►oco 

mox 

•‘XOO 

XO 

•00 

xoo 

CJ'  ►x 

X 

<iajo 

'JJXr^ 

ox  ► 

cro 

►LDX 

oox 

►  ►  I 

xxo 

►CJCU 

XXCM 

00  •• 

XO<l. 

O'  ••x 

0  *0 

X%T 

•‘XO 

<tx  •• 

0  -X 

XO'T 

oox 

►XX 

LMnJ-  ► 

►  Nj*r-4 

XO'CO 

ni  ••x 

X  ->r 

OQ  •- 

►XQ 

UJX  •• 

>rox 

X.OO 

mox 

►  ►^ 

l—lp^ 

•<30 

XXP** 

lU  ••x 

X 

n>in  •* 

OXLJ 

ox  •• 

ooox 

••CMX*»> 

oox 

•ox 

x<o— 

►^co  ► 

X  ••00 

xo  •‘X 

xoor*-^ 

xxo  •* 

CD  ►XO 

0  ••  ^X 

xo  •" 

OXO'T 

^  ^ 

oxx^ 

o^xx 

►Osr  - 

CL  X 

00-«'4“ 

►  ••x< 

x-r  -X 

••s: 

•  m  9^ 

► 

XXrOrO 

0  •‘XX 

xxo 

•» 

X  ••  •^JL 

ox  ••  •• 

mxcnx 

XI-  I-2.C 

<100 

UJ 

QXXO 

>rcDXX 

►r-CO 

x**^ 

h-4  •  ax^: 

U/LJU 

»- 

mooo 

••  ►XCJ 

XOOO  - 

.►a— 

0  ••  •>>—  •• 

<a: 

> 

••oo 

XXCMUJ 

irv  ►  ►x 

OfOlU 

«— •  —•  •x 

2:0 

cu 

X  ••  ••x 

OOO  ••o 

oxxx 

O  1  rM 

►2:  ^  ►i 

oxxo 

COQX  •• 

OCOXX 

► 

— W  wixx 

<i<t 

00000 

•OU'X 

►COOO 

-rcMKi 

X 

X 

►O' coo 

X  ••xr^ 

X  *0  •- 

X'-'l— 

**-iO  • 

h*  •* 

X 

XXOOO 

XX  -X 

x<x 

QCJ  (/>  •  ► 

►UJ 

oxxx 

xo  ••  •* 

xxxo 

X 

QC  •  1  •  •^ 

0 

OLUI- 

COiLG^ 

0-1-XX 

ocoo^r 

Xh- 

0  ••  •xz 

UJ<I 

uuk;>- 

•ox 

•*'•4*  ** 

XX'-' 

•  X-^XX— ' 

a:x 

00>tI3 

•*  ►  ► 

x^rm 

X  --x 

Ql— QC 

►•-•CMXXCO 

•- 

<aj 

<IXXX 

<i.rn  ••  •• 

<lOXCQ 

r<)'-'X 

<3L0^Xh-0 

o<i 

►GO 

»-ooo 

h-XXX 

^ujcacD 

A 

X  •Xh-  *0 

oixoi  — m 

<000 

<0^0 

^  1 

H-'-'X  ♦ 

01— r- <101'-^ 

LO  ►  ►  ► 

LJ  ••Xf-l 

UJ  ►  <•  ► 

IM  -tt  **■ 

w  UJ^  w 

>-f - CL 

XjuX 

X  ••  •* 

XXX 

‘—00  ••o 

1— 1-  I--  ••* 

CL  < 

CDUJULX 

cjoooo 

c/>r-xx 

zooxo 

<x  II  rno  II 

GQo  If  aooQO 

CLUJ 

Xi-^tU 

QCUJXCO 

ox»^ao 

hhCJOX 

HX  II 

00— U30X-2 

CL^h- LJh- 

0000 

LJOLnOO 

COOOO 

ujrsi^x 

VJXCMCXC/IO 

••<5L 

X 

W  u_ 

1—  UJ 

UUJ 

U.LJ 


it  -Hr  ’ft'  -ft-  ^ 


O 

U 


LJ 

dL 

UJ 


97 


•» 

-  X 

•» 

•k 

•kX 

X  ••o 

••x 

XX  •» 

xxu? 

xxo 

•►xcux 

0300 

OO^'T  •* 

xoox 

•-00  ^ 

-X 

••x 

XO  •Mn 

X  •ox 

•"XXO 

•k 

dj  »xo 

Uix  **0 

xu-o^ 

xxo^r 

oxor*  •k 

□JOXO 

UJ(Y)^  ^ 

sOOX  •* 

••xcMx 

ooo  •* 

(T^  ^  ••x 

UIOJOX 

xm  •‘04 

••  •ox 

••xxa3 

OO  ••x 

•‘xx 

XX  *<30 

xcMor- 

•k  .^XCD 

r^-X-^-O 

OOXUI 

r^u'< 

xxog  •• 

•kgr\<r  Vk 

ooooo 

CL)  •^x 

tncoxx 

XLOOX 

••o<i  ^ 

ox 

XL300 

GO  •»  *<U 

X  -ox 

•^ixtu 

oo 

V'XXX 

ox  -o 

xo>ro 

•*  •kX  •• 

•►xx  •* 

<10.0.0 

UJOO'  •• 

xxox 

xxxx 

00300 

Q>  ^  ^x. 

OO'TO 

XOOX 

•oo 

•-XXO 

k^sj-  •k^ 

m  wk  •kO 

X  ox 

XXO'^* 

•k  •kX  •* 

oxxo 

ax  o 

ar»>>r  •• 

XXsOX 

•“TX 

COOXO 

CO  •*  ••x 

XOCVJX 

OOO 

•‘XXX 

ox  ^oo 

xo 

^  •ox 

xmcuvu 

00X0 

<  •.X'L3 

J.J.  *<o 

r^oorsi  •• 

•k  •ktj  •. 

ox>u  •• 

00X0 

r*.  m>x 

xxxx 

•kO'^X 

OOUGOO 

^XXni 

0>T0«' 

XXOOl  •• 

••OCT^  •»  •» 

Xa'C7>CD 

XX  ••U' 

•>  9*^ 

—  X 

X  •.  .-X-- 

x<coo 

oox  - 

XXX  •‘X 

XX 

OXXOX 

coo  ^ 

-  -ox  •* 

QXXXX 

XX**'*^ 

0X00^ 

o  -^xxx 

xx^o-^ 

•MX  XXX 

2r 

OCDOOO 

•-XXOO 

oo  •->rx 

XOOCDO 

o 

•oo 

xxrn>r>r 

>4“-rx  •MX 

CD  •<>  •* 

—  •  XX 

X  -"X 

£X3(0 

-  •‘r-xir> 

XXX  •‘X 

X  •4JJh- 

>- 

•  •k 

oxxox 

(Tj  ••X  XX 

XX<IU'  •• 

Orn*^XCM 

•-^  •k^l— 

u 

•* 

CDOOOO 

•-xr-ooo 

oooxx 

•Mnr-xo 

o— ‘Z  d  •k 

icm 

<tco 

OOO  ••o 

xoooo 

OO)  •‘OO 

X  ^  * 

•Zw  Q. 

DO 

Xoo 

•»  •*  ^x  ^ 

CD  vQ  **^3  •• 

O  ••x 

OXX  •‘X 

wUJ  •kX 

u.  • 

f-UU 

xxxoox 

lU  •‘X  ••x 

-xr-x  •• 

xr-inxx 

•kUJxxx 

t-H 

•» 

oocx)o>r 

ox<xa 

XCJOXX 

•HX<<LX 

ox 

<XCJ 

uu 

OVP 00003 

•H^^XCX'O 

XXOV/'OO 

XOOU'O 

CJO 

Z^oo  • 

<v 

XQ 

h- 

X00O<O 

xo  •*  •‘X 

X  •»  •»  •*  •• 

fO  •• 

wQitO  d 

»- 

w< 

>-  •• 

XXX  XX 

X  •o  •• 

o  •‘XXO 

XXXXX 

z 

O  •  0  •k 

<x 

CD  ••UJ 

OUJXOO 

XXX  •‘X 

--'XOO 

oxxxx 

X  1  X 

0<  •k— 

UJ— 

UJH- 

OOoQoQOO 

oxxxo 

o^^•^^-x 

‘"'XXXX 

X  ivl 

O  •‘XZ 

K  tl 

a:< 

<►-> 

wOOOO 

wcnrooNj- 

•*  •<> 

oooo 

|f— .-"k  •* 

•  xxx-' 

zv 

DX 

X>CD 

•»  •» 

h»  ••xx ‘jr 

^  ^  •k  •k  •k 

OJ>^ 

H-iXXO 

C-Jh- 

CU 

<IXXXX 

<IXX  ••x 

<X(MO  •* 

)-XXXX 

Xww 

•kOXh"<C 

oo 

tu  •» 

h-OOOO 

h»r^xxo 

LJXXOX 

<t(\ICMXX 

o;kx 

<i  •:> 

Q.  • 

o^ujuj— **rn 

<oooo 

ca30<tO 

'TOOO 

O0‘^''XX 

m^x 

•W  • 

o 

V^X^-V)0|-^ 

^ 

l_»  •*  ^00  •• 

QJ  •»  Vk  9^^ 

•k  •‘OO 

1  t/> 

o 

CCw, 

>-<wCL 

xxxx 

XX  ••x 

xxxx  •• 

xxx  ••  •• 

/\X<fh  •••• 

I-  h- 

z 

tx  1 

□□ciju.2: 

0000300 

oooxin 

QsT— •oox 

l-O-MXX 

11  moo 

ooooo 

»<-4 

LU 

H-ILU 

Q.JJXXU. 

<iQXXO 

XXQ^^'^- 

<o<txrn 

'^UJ  tl  z  il 

ODQOO 

>— 

ooooo 

->00V'0 

cyi  ooooo 

QCooonoo 

<IX^X<L 

<1. 

LJ 

4^-J 

^  LU 

— 1— l-J-JO 

-J 

— 

U) 

xoooo 

u. 

U.UJ 

X  o 

<A<1<L<L<1 

»-«  X 

ooooo 

'WIT  ^ 


98 


ENC  n/MFUNC 


FLC^TING  POINT  TRIGONOMETRY  ROUTINE 
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ELSE  CALL  ACD  ( .PIONEHF  i  THA»THA  )  ; 
RETURN;  END; 


<X 

X 


•  «» 

•• 

< 

<l 

X 

X 

h- 

r-— 

«» 

X 

X 

XH 

oi 

1— 

^<l 

o<i 

♦-•X 

•-•X 

QlK 

a.K 

• 

•  •» 

-J 

X 

CD<t 

cm— 

ZJJ~ 

3  •" 

t/>*-4 

</)»-♦ 

0. 

•ft*  X 

O-W*  —1  # 

3  • 

0  -J^ 

1-4^ 

>-<OQ  •• 

lU  QQ 

OQ 

Xl-2<2: 

l-Z<Z 

J-ZX  UJ 

ZUJ  X 

<LX-J 

<1X3 

xcri--j 

Ctl-3 

oa  <t  ••* 

a  <1  ••* 

ao<t(\JoZ 

<iouz 

II  GL 

OL 

VC^*-hUJ«J 

o 

(/)^• 

••  coH- 

2: 

O-R-  LI _ \UJ 

LU 

OJ'VI-^UJU. 

LJ  ^*-4UJX 

tAJ 

U1 

X 

3 

>-< 

IXJ 

o 

tu 


102 


NONl:  PPCCECLPE  (PUNC,INFC); 
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IF  ((>(2)  AND  7FH)>50F)  THEN  GC  TO  STAR; 
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